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We propose a gauge-invariant and manifestly UV finite resummation of the 
physics of hard thermal/dense loops (HTL/HDL) in the thermodynamics of 
the quark-gluon plasma. The starting point is a simple, effectively one-loop 
expression for the entropy or the quark density which is derived from the fully 
self-consistent two-loop skeleton approximation to the free energy, but sub- 
^ ■ ject to further approximations, whose quality is tested in a scalar toy model. 

In contrast to the direct HTL/HDL-resummation of the one-loop free energy, 
in our approach both the leading-order (LO) and the next-to-leading order 
(NLO) effects of interactions are correctly reproduced and arise from kine- 
matical regimes where the HTL/HDL are justifiable approximations. The 
LO effects are entirely due to the (asymptotic) thermal masses of the hard 
particles. The NLO ones receive contributions both from soft excitations, as 
described by the HTL/HDL propagators, and from corrections to the dis- 
persion relation of the hard excitations, as given by HTL/HDL perturbation 
theory. The numerical evaluations of our final expressions show very good 
agreement with lattice data for zero-density QCD, for temperatures above 
twice the transition temperature. 
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I. INTRODUCTION 



Besides its obvious relevance for cosmology, astrophysics or ultra-relativistic heavy ion 
collisions, the study of QCD at high temperature and/or large baryonic density [|l],0 presents 
exciting theoretical challenges. It offers opportunity to explore the properties of matter in 
a regime where, unlike in ordinary hadronic matter, the fundamental fields of QCD — the 
quarks and gluons — are the dominant degrees of freedom and the fundamental symmetries 
are explicit. 

Unfortunately, analytical tools available for such a study are not many. However, because 
of asymptotic freedom, the gauge coupling becomes weak at high temperature, which invites 
us to try a perturbative treatment of the interactions. But explicit perturbative calculations 
of the QCD free energy at high temperature, which have been pushed in recent years up to 
the order ct'J 2 show an extremely poor convergence except for coupling constants as 

low asa < 0.05, which would correspond to temperatures as high as > 10 5 T C . Already the 
next-to-leading order perturbative correction, the so-called plasmon effect which is of order 
o?J 2 (x g 3 , signals the inadequacy of the conventional thermal perturbation theory except for 
very small coupling, because in contrast to the leading-order terms it leads to a free energy 
in excess of the ideal-gas value. 

Lattice results on the other hand show a slow approach of the ideal-gas result from 
below with deviations of not more than some 10-15% for temperatures a few times the 
deconfinement temperature. Besides, these results can be accounted for reasonably well by 
phenomenological fits involving massive "quasiparticles" [|5]|| with masses of the order of 
the perturbative leading-order thermal masses. This suggests that the failure of ordinary 
perturbation theory may not be directly related to the non-perturbative phenomena expected 
at the scale g 2 T and which cause a breakdown of the loop expansion at order g 6 and higher 
|jlj]. Rather, the quasiparticle fits support the idea that one should be able to give an 
accurate description of the thermodynamics of the QCD plasma in terms of its (relatively 
weakly interacting) quasiparticle excitations. 

It is worth emphasizing at this stage that, among the relevant degrees of freedom, the 
soft collective ones, with momenta of order gT, are already non-perturbative. Although 
their leading order contribution oc g 3 to the pressure can be easily isolated flT], it does not 
make much physical sense to regard this contribution as a genuine perturbative correction. 

Indeed, to leading order in g, the dynamics of the soft modes is described by an effective 
theory which includes the one-loop thermal fluctuations of the "hard" modes with momenta 
~ T. The relevant generalization of the Yang-Mills equation reads : 

D V F^ = ihlJ— — E l = W*J% + - + ... (1.1) 

where the induced current in the right hand side describes the polarization of the hard 
particles by the soft colour fields in an eikonal approximation. [In this equation, rho ~ gT 
is the Debye mass, E % a is the soft electric field, = (1, v), and the angular integral / dQ runs 
over the orientations of the unit vector v.] This current is non-local and gauge symmetry, 
which forces the presence of the covariant derivative = + igA^ in the denominator 



of Eq. (1.1), makes it also non-linear. When expanded in powers of A^, it generates an 



infinite series of non-local self-energy and vertex corrections, known as "hard thermal loops" 
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(HTL) 0,0). The latter encompass important physical phenomena, like screening effects and 
non-trivial dispersion relations for the soft excitations (and references therein). Similar 
phenomena exist also in the case of soft fermions, which, to leading order in g, obey the 
following generalized Dirac equation [p]] (with M ~ gT and ff = ■y fl v^) : 

^ = ^I^i^D)^ " ^ + W + - (1-2) 

At soft momenta k < gT, all HTL's are leading order effects, as obvious in Eqs. ( |1 . 1|) 
and ( |1.2|) , and must be consistently resummed. Analogs of HTL's exist at finite chemical 
potential /i. In the regime p ^> T these are often referred to as "hard dense loops" (HDL). 

In traditional perturbative calculations of the thermodynamics performed in imaginary 
time H, the HTL's play almost no role: only the Debye mass m 2 D needs to be resummed in 
the static electric gluon propagator [Kj. This resummation is responsible for the occurrence 



of odd powers of g in the perturbative expansion. 

Such a simple resummation however may become insufficient whenever a more complete 
information on the quasiparticles needs to be taken into account. Quite generally, this 
physical information is contained in the spectral weight p(k , k) related to the corresponding 
propagator by: 

D ^k)= r^fp^R. (1 . 3) 

j-oo 2tt k — to 

In the imaginary time formalism, and for bosonic fields, to = iuo n = i2irnT with integer n. 
Clearly, the restriction to the Matsubara mode with n = retains in the propagator only 
one moment of the spectral weight. In the HTL approximation, we know that the spectral 
density is divided into a pole at time-like momenta and a continuum at space-like momenta. 
While there exist physical observables which can be accurately described in perturbation 
theory by a single moment of the spectral weight, this does not appear to be the case in the 
calculations that we shall present and in which the various pieces of the spectral functions 
contribute in different ways. 

In fact, since the thermodynamical functions are dominated by hard degrees of freedom, 
an important effect of the soft modes will be to induce corrections on the hard quasiparticle 
dispersion relations. As we shall find, the spectral functions for large momenta will take 
the approximate form p(u, k) m 5{uj 2 — k 2 — m^J, where m 2 ^ ~ g 2 T 2 is the leading-order 
thermal mass (or asymptotic mass) of the hard excitation. Clearly, such an effect does not 
naturally emerge in a scheme where one resums just the n = Matsubara mode. 

In order to overcome all these limitations, it has been recently proposed to perform full 
resummations of the HTL self-energies U^ u and £ in calculations of the thermodynamical 



functions. In Refs. |TT|,[T2[, this has been done by merely replacing the free propagators by 
the corresponding HTL-resummed ones in the expression of the free-energy of the ideal gas; 
e.g. (in simplified notations) : 

Tr log £>„ 1 — > Tr log^ 1 + n) . (1.4) 

In principle, this is just the first step in a systematic procedure which consists in resumming 
the HTL's by adding and subtracting them to the tree-level QCD Lagrangian. This would 
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be the extension to QCD of the so-called "screened perturbation theory" [13, 14 1 , a method 



which, for scalar field theories, has shown an improved convergence (in one- and two-loop 
calculations) as compared to the straightforward perturbative expansion. But in its zeroth 
order approximation in Eq. (|1.4|) , this method over-includes the leading-order interaction 
term oc g 2 (while correctly reproducing the order-g 3 contribution), and gives rise to new, ul- 
timately temperature-dependent UV divergences and associated additional renormalization 
scheme dependences. 

Another drawback of such a direct HTL resummation appears to be that the HTL's are 
kept in the hard momentum regime where they are no longer describing actual physics, while 
hard momenta are providing the dominant contributions to the thermodynamic potential. 

Our approach on the other hand [^,|l^] will be based on self-consistent approximations 
using the skeleton representation of the thermodynamic potential JT7J which takes care 
of overcounting problems automatically, without the need for thermal counterterms. We 
shall mainly consider the so-called 2-loop-$-derivable ]TB| approximation, for which it turns 
out that the first derivatives of the thermodynamic potential, the entropy and the quark 
densities, take a rather simple, effectively one- loop form |19|p0| , but in terms of fully dressed 
propagators. 

In gauge theories, the generalized gap equations that determine these dressed propaga- 
tors are too complicated to be solved exactly (even numerically). But an exact solution 
would anyhow be unsatisfactory because <I>-derivable approximations in general do not re- 
spect gauge invariance. We therefore propose gauge independent but only approximately 
self-consistent dressed propagators as obtained from (HTL) perturbation theory. Using 
these in the entropyf] expression obtained from the 2-loop-$-derivable approximation gives 
a gauge-independent and UV finite approximation for the entropy, which, while being non- 
perturbative in the coupling, contains the correct leading-order (LO) and the next-to-leading 
order (NLO) effects of interactions in accordance with thermal perturbation theory. Both 
turn out to arise from kinematical regimes where the HTL's are justifiable approximations. 

While also being effectively a resummed one-loop expression, the approximately self- 
consistent entropy differs from the direct HTL-resummation of the free energy in Eq. ( |1.4|) 
in that it includes correctly also the LO interaction effects. Remarkably, in our approach the 
latter are entirely determined by the (asymptotic) thermal masses of the hard excitations. 
This agrees with and justifies the simple quasiparticle models of Ref. HH, which assume 
constant masses equal to the respective asymptotic thermal masses for quarks and as many 
(scalar) bosons as there are transverse gluons. Whereas these models do not include the 
correct NLO (plasmon) effect, our approach does, but in a rather unconventional manner 
which demonstrates the nontriviality of the resummation that has been achieved: only part 
of the plasmon effect is coming directly from soft excitations; a larger part arises from 
corrections to the dispersion relation of the (dominant) hard excitations by soft modes, as 
determined by standard HTL perturbation theory H . 

Because of the approximations that we have made, it does matter whether the entropy 
or the thermodynamic potential is considered. Our approach however attempts to take 



For brevity we refer only to the entropy explicitly, but all of the following remarks apply to the 
density as well. 
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advantage of the fact that entropy is generally the simpler quantity. Indeed, the way by 
which the LO and NLO interaction contributions can be traced to spectral properties of free 
quasiparticles within our entropy expressions indicates a posteriori the adequateness of this 
particular resummation scheme to the physics contained in the HTL propagators. 

The present paper is organized as follows: In Sect. II, the general formalism of ^-derivable 
self-consistent approximations is reviewed and the central, effectively one-loop formula for 
the entropy in a two-loop skeleton approximation to the thermodynamic potential is derived 
in a scalar theory with cubic and quartic interactions. In the simple solvable model of 
large- N scalar O(N) theory |ZI,2"!Z|], where the two-loop ^-derivable approximation becomes 



exact, the further approximations that will be considered in the QCD case are compared 
with the exact solution and their renormalization scale dependence is exhibited. 

In Sect. Ill, the approximately self-consistent resummations are introduced for purely 
gluonic QCD first, and equivalence with conventional perturbation theory up to and in- 
cluding order g 3 is proved and analyzed in detail. Sect. IV generalizes this to QCD with 
quarks and to the quark density as an additional thermodynamic quantity. Some of the 
more technical details of how the plasmon effect arises in our approach are relegated to the 
Appendix. 

In Sect. V, the various approximations are evaluated numerically. We find that the 
plasmon effect, which is largely responsible for the poor convergence properties of conven- 
tional thermal perturbation theory, in our approach leads only to moderate contributions 
when compared with the leading-order effects. When combined with a two-loop renormal- 
ization group improvement, our results are found to compare remarkably well with available 
lattice data for temperatures above twice the deconfinement temperature. Moreover, we 
also present numerical results for the quark density at zero temperature and large chemical 
potential. 



II. GENERAL FORMALISM. THE SCALAR FIELD 

In this section we develop the formalism of propagator renormalization using tech- 
niques that allow systematic rearrangements of the perturbative expansion avoiding double- 
countings. We shall recall in particular how self-consistent approximations can be used to 
obtain a simple expression for the entropy which isolates the contribution of the elementary 
excitations as a leading contribution. To get familiarity with the formalism, we demonstrate 
some of its important features with the example of the scalar field. This provides, in par- 
ticular, a test of the validity of approximations which will be used in dealing with QCD in 
the rest of the paper. 



A. Skeleton expansion for thermodynamical potential and entropy 



The thermodynamic potential Q = —PV of the scalar field can be written as the following 
functional of the full propagator D ||i7| , |i"8| : 



PQ[D] 



log Z = - Tr log D' 1 - i Tr IXD + $[£>] 



(2.1) 
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where Tr denotes the trace in configuration space, (3 = 1/T, II is the self-energy related to 
D by Dyson's equation (D denotes the bare propagator): 



zr 1 = d x + n, (2.2) 

and $[-D] is the sum of the 2-particle- irreducible "skeleton" diagrams 

-$[£>]= 1/12^^ +1/8 QQ +1/48^^ +... (2.3) 

The essential property of the functional Q[D] is to be stationary under variations of D 
(at fixed Dq) around the physical propagator. The physical pressure is then obtained as the 
value of Q[D] at its extremum. The stationarity condition, 

5Q[D]/5D = 0, (2.4) 

implies the following relation 

6$[D]/6D = -II, (2.5) 

which, together with Eq. ( j2.2|) , defines the physical propagator and self-energy in a self- 
consistent way. Eq. (|2.5|) expresses the fact that the skeleton diagrams contributing to II are 
obtained by opening up one line of a two-particle-irreducible skeleton. Note that while the 
diagrams of the bare perturbation theory, i.e., those involving bare propagators, are counted 
once and only once in the expression of II given above, the diagrams of bare perturbation 
theory contributing to the thermodynamic potential are counted several times in $. The 
extra terms in Eq. (\2.1\ ) precisely correct for this double-counting. 

Self-consistent (or variational) approximations, i.e., approximations which preserve the 
stationarity property Q2.4j) , are obtained by selecting a class of skeletons in $[-D] and calcu- 
lating II from Eq. ( |2.5| ). Such approximations are commonly called "<I>-derivable" fTSf . 

The traces over configuration space in Eq. fl2.1|) involve integration over imaginary time 
and over spatial coordinates. Alternatively, these can be turned into summations over 
Matsubara frequencies and integrations over spatial momenta: 

J dr J d 3 x -> /3V J [dk], (2.6) 

where V is the spatial volume, k^ = (iu n , k) and u n = rnrT, with n even (odd) for bosonic 
(fermionic) fields (the fermions will be discussed later). We have introduced a condensed no- 
tation for the the measure of the loop integrals (i.e., the sum over the Matsubara frequencies 
uj n and the integral over the spatial momentum k): 

Strictly speaking, the sum-integrals in equations like Eq. ( |2.1| ) contain ultraviolet diver- 
gences, which requires regularization (e.g., by dimensional continuation). Since, however, 
most of the forthcoming calculations will be free of ultraviolet problems (for the reasons 
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explained at the end of this subsection), we do not need to specify here the UV regulator 
(see however Sect. II B for explicit calculations). 

For the purpose of developing approximations for the entropy it is convenient to perform 
the summations over the Matsubara frequencies. One obtains then integrals over real fre- 
quencies involving discontinuities of propagators or self-energies which have a direct physical 
significance. Using standard contour integration techniques, one gets: 

r d 4 k 

Q/V = J-^—^n(uo)(\m\og(-uj 2 + k 2 + Yi)- ImUDj + T<^[D]/V (2.8) 

where n(uo) = l/(e^ - 1). 

The analytic propagator D(u,k) can be expressed in terms of the spectral function: 

D^h)= rpihlV. (2.9) 

J-oo 27T ko — U 

and we define, for u> real, 

lmD(u,k) = Im D(u + ie,k) = p ^ k \ (2.10) 

The imaginary parts of other quantities are defined similarly. 
We are now in the position to calculate the entropy density: 

S = -d(n/V)/dT. (2.11) 



The thermodynamic potential, as given by Eq. ( j2.8| ) depends on the temperature through 
the statistical factors n(u) and the spectral function p, which is determined entirely by the 
self-energy. Because of Eq. (|2.4j) the temperature derivative of the spectral density in the 
dressed propagator cancels out in the entropy density and one obtains [^9j , p0[| : 

s = - f?j dn } u) ImlogZrV,fc) 
J (2tt) 4 dT B v ' ; 

+ /(0i^ Imn (^ fc ) Re ^(^ fc ) + 5, ( 2 - 12 ) 



with 



S' 



d(T$) r d 4 k dn{uj 



dT 

We shall verify explicitly that for the two-loop skeletons, we have: 

S' = 0. (2.14) 

Loosely speaking, the first two terms in Eq. ( |2.12[ ) represent essentially the entropy of "in- 
dependent quasiparticles" , while S' accounts for a residual interaction among these quasi- 
particles [20| . 
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Since the condition (|2.14|) plays an important role in our work, we shall derive it explicitly 
in a scalar model with interaction term 

C mt = {g/m* ~ (A/4!)0 4 , 

which is a simple toy model of the tri- and quadrilinear self-interactions of gauge bosons. 
(Interactions with fermions are already covered by the analysis contained in Ref. PD| .) In 
the two-loop approximation, where only the first two diagrams of the skeletons in Eq. ( f2.3|) 
are kept, the contribution involving two 3- vertices reads 



12 



T 2 E 



wi,a>2 



d 3 ki d 3 k 2 
(2tt) 6 



Dfa, \h\)D{u 2 , IhW-U! -co 2 ,\-h- h\). (2.15) 



Expressing the propagators in terms of the spectral functions, and evaluating the Matsubara 
sums by contour integration, one gets: 



y<j>(») 



g_ 

12 



2 r d 4 kd 4 k'd 4 k" 



■5 3 (k + k'+ k")p(k)p(k')p(k")P- 



l 



(27T) 9 ' V ~ ' ' ^vvrv/rv /- ko + k/Q + k ,, 

x {[n(ko) + l][n%) + n(A#) + 1] + n{k' )n{k^)} 



where P denotes the principal value prescription and we have used the identity: 

n(x +y)[l + n(x) + n(y)] = n(x)n(y) . 
The two-loop skeleton involving the 4-vertex is given by the simpler expression 



(2.16) 



(2.17) 



_ ycjj(f>) 



A 



E 



d 3 k 



A rd 4 kd A k' 



(27T) 



p(k)p(k'){n(k )n(k> )}. (2.18) 



According to Eq. ([2.13|) , the first contribution to S' is given by differentiating Eqs. (p. 16 ) 
and ( J2.18 ) with respect to T at fixed p. Because the integrand in front of the curly brackets 
in (|2.16|) is symmetric, the arguments of the distribution functions can be freely exchanged 
as long as the fact that their products come with distinct arguments is preserved. dx[— T$] 
is therefore obtained by replacing the terms in curly brackets in ( [2.16| ) by {6n(k' )dTn(k ) + 
3d T n(k )}} and that in (glg|) by 2n{k' )d T n{k ). 

The second contribution to S' involves the real part of the self-energy as given by the two 
(dressed) one- loop diagrams following from opening up one line in the first two diagrams in 

£3D, 



Re n<»V, 9 ) - 4 \i + ,1) 



x[n(k ) + n(k' ) + 1]P 



to + k + k' 



ReUV = ±J-^n(k )p(k ,k) 



(2.19) 
(2.20) 



This gives 
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d'k ^0) Ren (a) ImD = _I [ d4kd * k ' f k " 6 3 { k + k> + k")p{k)p{k') P {k") 

4 J (27T) y 



J (2tt) 4 dT A J (2tt 

* P ko + V +k A d Mko)}ln(k' )+n(k>>) + 1] (2.21) 

(2 - 22) 

where we have used ImD = p/2. Indeed, this cancels precisely — <9t[T$] as obtained above, 
verifying the proposition that 5' = for the lowest-order (two- loop) diagrams in $[-D]. 

As the previous derivation shows, the vanishing of S' holds whether the propagator are 
the self-consistent propagators or not. That is, only the relation ( |2.5| ) is used, and the proof 
does not require D to satisfy the self-consistent Dyson equation ( p.2|) . A general analysis of 
the contributions to S' and their physical interpretation can be found in Ref. |23| . 



We emphasize now a few attractive features of Eq. Q2.12j ) with <S' = 0, which makes the 
entropy a privileged quantity to study the thermodynamics of ultrarelativistic plasmas. We 
note first that the formula for S at 2-loop order involves the self-energy only at 1-loop order. 
Besides this important simplification, this formula for S, in contrast to the pressure, has the 
advantage of manifest ultra-violet finiteness, since dn/dT vanishes exponentially for both 
uj — > ±oo. Also, any multiplicative renormalization D — > ZD, U — > Z~ l Tl with real Z drops 
out from Eq. ( |2.12| ). Finally, the entropy has a more direct quasiparticle interpretation than 
the pressure. This will be illustrated explicitly in the simple model of the next subsection. 
More generally, Eq. ( |2.12| ) can be transformed with the help of the following identity: 

/ j m TT \ 

Im log IT 1 (u, k) = arctan - — - ne(u)6{- Re D' 1 ), (2.23) 



ReD 



with e(u) the sign function and — | < arctan(x) < |. Using this identity we rewrite S as 

<S = Sp l e + Sdampi with 

d 3 k 

(2^)3 {(! +n k )\og(l + n k ) - n k \ogn k ). (2.24) 

To get the second line, we have made an integration by part, using 
dn(uo) dcr(uj) 



dT duj 



a(u) = -n\ogn + (1 +n)log(l+n), (2.25) 



and we have set n k = n(e k ), with e k solution of ReZ) -1 ^ = e k , k) = 0. The quasiparticles 
thus defined by the poles of the propagator are sometimes called "dynamical quasiparticles" 
pSfl . The quantity S po i e is the entropy of a system of such non-interacting quasiparticles, 
while the quantity 



damp 



d 4 k dn{uj) 
(2tt) 4 dT 



, ImU 

Imil(u;, k) Re-D(a;, k) — arctan 



ReD" 



^2.26) 



which vanishes when Imll vanishes, is a contribution coming from the continuum part of 
the quasiparticle spectral weights. 
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B. A simple model 



In this section we shall present the self-consistent solution for the (A/4!)</> 4 theory, keeping 
in $ only the two- loop skeleton whose explicit expression is given in Eq. (|2.18| ). Anticipating 
the fact that the fully dressed propagator will be that of a massive particle, we write the 
spectral function as p(fcn,k) = 2vt e(ko) S(k 2 — k 2 — m 2 ), and consider m as a variational 
parameter. The thermodynamic potential (|2.1|) , or equivalently the pressure, becomes then 
a simple function of m. By Dyson's equation, the self-energy is simply II = m 2 . We set: 

/M^/[W)4/[d ^ + J 2 + m , (2.27) 

Then the pressure can be written as: 

„ il lr d 3 k 1 r d 3 k , fl _ . 9t/ , An t9/ n ,„ „ M 

- P = 7= 2/ (2?F £i + ^/ (2^ l0g(1 " 6 » " m I(m) + i ' (m) ' (2 - 28) 

where e\ = k 2 + m 2 . By demanding that P be stationary with respect to m one obtains the 
self-consistency condition which takes here the form of a "gap equation": 

m 2 = A /(m). (2.29) 

The pressure in the two-loop ^-derivable approximation, as given by Eqs. fl2.27| )-( f2.29|) , is 
formally the same as the pressure per scalar degree of freedom in the (massless) A-component 
model with the interaction term written as (A/4!)(0j0,) 2 in the limit A — > oo [22]. From 



the experience with this latter model, we know that Eqs. (|2.27| )- (|2.29|) admit an exact, 
renormalizable solution which we recall now. 

At this stage, we need to specify some properties of the loop integral J(m) which we can 
write as the sum of a vacuum piece Io(m) and a finite temperature piece Iri^n) such that, at 
fixed m, Ixi^n) —> as T — > 0. We use dimensional regularization to control the ultraviolet 
divergences present in 1$, which implies /o(0) = 0. Explicitly one has: 

///(m) = ( - + log 4 + 1 )+ Mm) + 0(e), (2.30) 



32vt 2 \ e m 



with 



f d 3 k n(e k ) 

/T(m)= i(27F^r' (2 - 31) 

ande fc = (k 2 + m 2 ) 1 / 2 . In Eq. ( g]§ ), \x is the scale of dimensional regularization, introduced, 
as usual, by rewriting the bare coupling An as //An, with dimensionless Ao; furthermore, 
e = 4 — n, with n the number of space-time dimensions, and ft 2 = 4vTe _7 /i 2 . 

We use the modified minimal subtraction scheme (MS) and define a dimensionless renor- 
malized coupling A by: 

111 

A ~ A^ + 16^' (2 ' 32) 
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When expressed in terms of the renormalized coupling, the gap equation becomes free of 
ultraviolet divergences. It reads: 



m 



A f d 3 k n(e k ) ^ Am 2 



(27T) 



32tt 5 



mr 



(2.33) 



The renormalized coupling constant satisfies 

dA A 2 



dlog [L 



16tt 2 



^2.34) 



which ensures that the solution m 2 of Eq. ( p. 33 ) is independent of p,. Eq. ( 2.34 ) coincides 
with the exact /5-function in the large- N limit, but gives only one third of the lowest-order 
perturbative /5-function for N — 1. This is no actual fault since the running of the coupling 
affects the thermodynamic potential only at order A 2 which is beyond the perturbative 
accuracy of the 2-loop ^-derivable approximation. In order to see the correct one-loop 
/5-function at finite N, the approximation for $ would have to be pushed to 3-loop order. 

Note also that, in the present approximation, the renormalization ( |2.32| ) of the coupling 
constant is sufficient to make the pressure ( |2.28|) finite. Indeed, in dimensional regularization 
the sum of the zero point energies £fc/2 in Eq. Q2.28p reads: 



Jn-li 



(2rr) n - 1 2 



m 



64vr 2 



- + log^ + |)+0(e) 
e m z 2 1 



(2.35) 



so that 



in— 1 



k e k 



(27I-)"- 1 2 



2A 



m 
2A 



m 

64^ 



log 



m 2 



3 



+ 0(6) 



(2.36) 



is indeed UV finite as n — > 4. After also using the gap equation ( |2.33| ), one obtains the 
/i-independent result 



P 



-T 



log(l 



1 4 

H I T (m) H -. 

2 v 1 128vr 2 



(2.37) 



We now compute the entropy according to Eq. fl2.12j ). Since Imll = and Re II = m 2 
we have simply: 



S 



d 4 k dn(u) 
(2tt) 4 dT 



Im log(A; 2 — ijj 2 + m 2 ). 



Using 



Im log(A; 2 - uo 2 + m 2 ) = -ire(uj)9(uj 2 - e\), 

and the identity ( |2.25| ), one can rewrite Eq. ( |2.38|) in the form (with n k = n(efc)): 

/" c/^ k c *i 
5 = y^p|(l + ri A; )log(l + n s .) - 7i fc lognfc]. 



(2.38) 



(2.39) 



(2.40) 
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This formula shows that, in the present approximation, the entropy of the interacting scalar 
gas is formally identical to the entropy of an ideal gas of massive bosons, with mass m. 

It is instructive to observe that such a simple interpretation does not hold for the pressure. 
The pressure of an ideal gas of massive bosons is given by: 

pm ^ - IS? f dw (" M + s) = - Iwf { Tlog(1 " ^' T) + f 1 • (2 ' 41) 

which differs indeed from Eq. ( p.28|) by the term m 4 / A which corrects for the double-counting 
of the interactions included in the thermal mass. Note that since the mass depends on the 
temperature, and since S = dP/dT, it is not surprising to find such a mismatch. 

Moreover, unlike the correct expression ( [2.2S| ), Eq. ( |2.41| ) is afflicted with UV divergences 
which in dimensional regularization are proportional to m 4 (cf. Eq. ( 2.35|) ), and hence 



dependent upon the temperature. This is precisely the kind of divergences which are met 



in the one-loop HTL-resummed calculation of the pressure in QCD of Ref. [I I 



C. Comparison with thermal perturbation theory 

In view of the subsequent application to QCD, where a fully self-consistent determi- 
nation of the gluonic self-energy seems prohibitively difficult, we shall be led to consider 
approximations to the gap equation. These will be constructed such that they reproduce 
(but eventually transcend) the perturbative results up to and including order A 3 / 2 or g 3 , 
which is the maximum perturbative accuracy allowed by the approximation S' = 0. 

In view of this it is important to understand the perturbative content of the self-consistent 
approximations for m 2 , P and S. In this section we shall demonstrate that, when expanded 
in powers of the coupling constant, these approximations reproduce the correct perturbative 
results up to order A 3 / 2 0. This will also elucidate how perturbation theory gets reorganized 
by the use of the skeleton representation together with the stationarity principle. 

For the scalar theory with only (A/4!)0 4 self-interactions, we write?] A = 24g 2 , and 
compute the corresponding self-energy II = m 2 by solving the gap equation ( |2.33|) in an 
expansion in powers of g, up to order g 3 . Since we anticipate m to be of order gT, we 
can ignore the second term oc Am 2 ~ g 4 in the r.h.s. of Eq. ( |2.33|) , and perform a high- 
temperature expansion of the integral Ir(w) in the first term (cf. Eq. ( |2.31| )) up to terms 
linear in m. This gives the following, approximate, gap equation: 

m 2 ~ g 2 T 2 - - g 2 Tm . (2.42) 

7T 

The first term in the r.h.s. arises as 

24/MO) = 12//^ ^ = g 2 T 2 = m 2 . (2.43) 



2 This normalization for g is chosen in view of the subsequent extension to QCD since it makes 
the scalar thermal mass in Eq. ( 2.43| ) equal to the leading-order Debye mass in pure-glue QCD 
(Eq. (|3l6|) with N = 3). 
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This is also the leading-order result for m 2 , commonly dubbed the "hard thermal loop" 
(HTL)| |||{]] because the loop integral in Eq. fl2.43|) is saturated by hard momenta k ~ T. 
The second term, linear in m, in Eq. ( |2.42| ) comes from 



12g 2 



d 3 k ( n{ek) n(k) 



(2*) 



~k 



Ylg 2 T I- 
J 



d 3 k 



(2tt) 3 VA; 2 + 



77T 



1 



3<T 



7T 



mT . 



(2.44) 



where we have used the fact that the momentum integral is saturated by soft momenta 
k ~ gT, so that to the order of interest n(ek) — T/e^ (and similarly n(k) ~ T/k). This 
provides the next-to-leading order (NLO) correction to the thermal mass 



5m 



3<r 



7T 



mT 



-9 3 T 2 . 



(2.45) 



Thus, to order g 3 , one has m 2 



in 



+ 5m 2 . In standard perturbation theory |],|2|], the 
first term arises as the one-loop tadpole diagram evaluated with a bare massless propagator, 
while the second term comes from the same diagram where the internal line is soft and 
dressed by the HTL, that is D(u, k) = — 1/ {uo 2 — k 2 — m 2 ). (At soft momenta k ~ m ~ gT, 
m 2 is of the same order as the free inverse propagator Dq 1 ~ k 2 ~ (? 2 T 2 , and thus cannot 
be expanded out of the HTL-dressed propagator D(u,k).) 

Consider similarly the perturbative estimates for the pressure and entropy, as obtained 
by evaluating Eqs. ( 2.28|) and ( 2.40| ) with the perturbative self-energy II = m 2 ~ m 2 + 5m 2 , 



and further expanding in powers of g, to order g . The renormalized version of Eq. (|2.28 ) 



yields, to this order (recall that m ~ gT and A 

7T 2 T 4 



9' 



P 



90 



m 2 T 2 m 3 T 

+ ~m + ■ 



24 



12tt 



4 

m 
+ 2V 



(2.46) 



The first terms before the dots represent the pressure of massive bosons, i.e. Eq. fl2.41p 
expanded up to third order in powers of m/T. From Eq. ( |2.46| ), it can be easily verified 
that the above perturbative solution for m 2 ensures the stationarity of P up to order g 3 , as 
it should. Indeed, if we denote 



m 2 T 2 
'~2T 



P 2 (m) 



+ 



m 
2A 



PM = 



then the following identities hold: 
dP 2 



dm 



0. 



OPo 



dm 



+ 



dP* 



rh+5m 



m 3 T 
12tt 



0. 



This shows that the NLO mass correction 5m 2 

(dP 3 /dm) 



5m 



(d 2 P2/dm 2 



dm 

g 3 T 2 can be also obtained as 
= m , 

7T 



(2.47) 



(2.48) 



^2.49) 



3 In the following, HTL quantities will be marked by a hat. 
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in agreement with Eq. fl2.45|) . Moreover, P2 = P2{m) = — g 2 T 2 /48 and P3 = P^{m) = 
m 3 T /12n are indeed the correct perturbative corrections to the pressure, to orders g 2 and 
3J. In fact, the pressure to this order can be written as: 



g 3 , respectively 



P 



vr^T 4 m 2 T 2 



90 

7T 2 r 4 

90 



24 
,2 



3 ^ 

7T 



m 3 T 
12tt 



th." „ 771 T 

— T + 

48 12tt 



(2.50) 



Note that the term of order g 2 is only half of that one would obtain from Eq. ( |2.41| ) by 
replacing m by rh. This is due to the aforementioned mismatch between Eq. (|2.41|) and 
the correct expression for the pressure, Eq. ( |2.28| ). In fact, going back to Eq. (2.1), one 
observes that the net order g 2 contribution to the pressure comes from $ evaluated with bare 
propagators: the order g 2 contributions in the other two terms mutually cancel indeed. This 
is to be expected: there is a single diagram of order g 2 \ this is a skeleton diagram, counted 
therefore once and only once in $. Observe also that the terms of order g 3 originating 
from the terms rh 2 and rh 4 mutually cancel; that is, the NLO mass correction 5m drops 
out from the pressure up to order g 3 . This is no accident: the cancellation results from the 
stationarity of P at order g 2 , the first equation 



Consider now the entropy density. The correct perturbative result up to order g may 



be obtained directly by taking the total derivative of the pressure, Eq. ( |2.50 ) with respect 
to T. One then obtains: 




m 2 T 2 m 3 T 



48 



+ 



12tt 



+ 0{g% 



(2.51) 



We wish, however, to proceed differently, using Eq. (|2.40|) , or equivalently, since 
dPj dm = when m is a solution of the gap equation, by writing: 



S 



OP 
df 



[2.52) 



This yields: 



„ 4 /tt 2 T 4 m 2 T 2 m 3 T\ 4 ^ 

S = + + 0(m 4 /T), 

T V 90 48 48tt ) V 11 



(2.53) 



which coincides as expected with the expression obtained by expanding the entropy of mas- 
sive bosons, Eq. ( |2.40| ), up to order (m/T) 3 . If we now replace m by its leading order value 
771, the resulting approximation for S reproduces the perturbative effect of order ~ g 2 , but it 
underestimates the correction of order g 3 by a factor of 4. This is corrected by changing m 
to rn. + 5m with 5m = —3gm/27r in the second order term of Eq. (|2.53 ). Note that although 
it makes no difference to enforce the gap equation to order g 3 in the pressure (because of 
the cancellation discussed above), there is no such cancellation in the entropy. 

In view of the forthcoming application to QCD, we shall now rephrase the previous 
discussion in slightly more general terms, though still restricted to the main simplification 
that the present simple model offers: a self-energy that is constant and real. 
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Because of the stationarity of the thermodynamic potential, Eq. ( |2.4| ), the order g 3 term 
in the pressure is coming entirely from the log term in the thermodynamic potential with 
II = ft, which reads: 



T 
~2 



d*k T 
(2vr) 4 w 
d 3 k 
(2vr) 3 



Im 



log(l + D (u, k)fl) - D (u, k)fl 



log 1 + 



n 

k 2 



n 

¥ 



12tt 



n 3/2 



(2.54) 



where we have subtracted the order-g 2 contribution and used the fact that the remain- 
ing integrand is dominated by soft momenta to replace n(uj) by T/u. The corresponding 
contribution to the entropy follows as: 



d^ 
dT 



dP* 



dT 



+ 



dU 



™ s5 

t dT 



3 



3 5 



(2.55) 



where , the derivative of P3 at constant IT, equals 1/4 of the total order-g 3 entropy. The 
remaining 3/4 come from the derivative of II. 

Alternatively, the entropy can be obtained from our master equation (|2.12|) which, in the 
present model where Im II = 0, simplifies into: 

d 4 k dn{uj 



S 



Im log-D 1 (u, k) 



(2?r) 4 dT 

The term of order g 2 is obtained by writing \0gD~ 1 = logr^" 1 +log(l + -D n), setting II 
and expanding the logarithm to first order in II. One then obtains: 

d 4 k dn(u>) 



2.56) 

n 



n Im D (u,k). 



(2.57) 



(2tt) 4 dT 

Since ImDo(w, k) = -Ke(u)5(u 2 — k 2 ), the integrand in (|2.57| ) is concentrated on the unper- 
turbed mass-shell. The ensuing momentum integral immediately yields S2 = —TIT/12, in 
agreement with Eq. (|2.51|) . 



:(b) 



According to Eq. ( [2.55] ), the contribution of order g 3 involves two pieces, S 3 = S3 +S : 
(cf. Eq. (|2.55| )). These can be also understood as the contributions to Eq. ( |2.56|) from 
different momentum regimes. Specifically, the soft momenta in the latter yield: 

d 4 k 1 



S: 



soft 



(2tt) 4 uo 



Im 



log(l + D U) - D U 



(2.58) 



which is the same as <Sg in Eq. fl2.55| ). The second contribution of order g 3 comes from 
hard momenta in Eq. ( |2.56| ), and is obtained by replacing II — > 5U in Eq. (|2.57Q . This yields 

d 3 k 1 dn(k) 



chard 
^3 



2 

t dn 

2" dT 



(2tt) 3 k dT 
d 3 k ( 1 



(2tt) 3 
d 4 k T 

Im 



(2vr) 4 uj 



\k 2 + U 
dfl 



A dT 
" k 2 ) 



dT 



(D-D 



^3 5 



(2.59) 



where we have used Eq. fl2.43|) for fl in the first line and Eq. (|2.45|) for SU in the second line. 
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D. Approximately self-consistent solutions 



As we have seen, the 2-loop ^-derivable approximation provides an expression for the 
entropy S as a functional of the self-energy II — namely, Eq. (|2.12| ) with S' = — which has 
a simple quasiparticle interpretation and is manifestly ultraviolet finite for any (finite) IT. 
These attractive features of Eq. Q2.12Q are independent of the specific form of the self-energy, 
and will be shown to hold in QCD as well. Of course, within this approximation, the self- 
energy is uniquely specified: by the stationarity principle, this is given by the self-consistent 
solution to the one-loop gap equation. In the scalar 4 -model, it was easy to give the exact 
solution to this equation (cf. Sect. 2.B), which coincides with the well-known solution of 
a scalar 0(iv)-model in the limit N — ► oo ||22|| . In QCD, however, it will turn out that a 
fully self-consistent solution is both prohibitively difficult (because of the non-locality of the 
gap equation), and not really desirable (for reasons to be discussed in Sect. 3.B below). 
This leads us to consider approximately self- consistent resummations, which are obtained 
in two steps: (a) An approximation is constructed for the solution II to the gap equation, 
and (b) the entropy (|2.12|) is evaluated exactly (i.e., numerically) with this approximate self- 
energy. While step (b) above is unambiguous and inherently nonperturbative, step (a), on 
the other hand, will be constrained primarily by the requirement of preserving the maximum 
possible perturbative accuracy, of order g 3 (cf. Sect. 2.C). In addition to that, we shall add 
the qualitative requirement that the approximation for IT, and the ensuing one for S, are 
well defined and physically meaningful for all the values of g of interest, and not only for 
small g — that is, for all the values of g where the fully self-consistent calculation makes 
sense a priori. As we shall shortly see, this last requirement generally excludes a strictly 
perturbative solution to the gap equation. 

Of course, even with this last requirement, there is still a large ambiguity in the choice of 
the approximate self-energy. In this respect the scalar 4 -model provides an opportunity for 
testing the quality of these approximations against the exact solution of the gap equation of 
the fully self-consistent two-loop calculation. Similar approximations will be subsequently 
used in QCD. 

The exact solutionfj of the gap equation is determined by the transcendental Eq. fl2.33| ) 
with A = 24g 2 . With p, = 2ttT, the result m/T as a function of g is given by the full line 
in Fig. [I]. As an exact result, it is independent of the renormalization scale: a change of 
/2 —>■ p,' has to be followed by a change of the renormalized coupling g(fl) —>■ g(pf) according 
to (cf. Eq. (pj))g 



4 More precisely, as discussed in detail in Ref. [22|, Eq. ( 2.33| ) has two solutions, a fact that is 



frequently overlooked. The larger of the two is exponentially larger than T for small coupling and 
has to be ruled out because our scalar model is consistent only as an effective (cut-off) theory. 

5 So the scalar theory is fully defined by giving both a dimensionful scale \x and the associated 
coupling strength g(fi). Equivalently, as usually done in QCD, we could just give a scale and 
agree e.g. that ff(A^) = +00. In this section we shall take the former point of view, so for any 
given temperature T, different values of g(27rT) parametrize differently coupled theories. 



17 



m/T 




FIG. 1. Comparison of the exact thermal mass in the large- ./V scalar 0(iV)-model as a function 
of g(p = 2ttT) (full line) with the leading-order (HTL) perturbative result (long-dashed lines), 
the NLO one corresponding to rh 2 + 5m 2 (shorter-dashed lines), and the perturbatively equivalent 
(NLA) version ( p. 62 ) (dotted lines). Except for the exact result, all these are renormalization scale 
dependent, the central lines corresponding to p! = 2irT, the adjacent ones to ttT and 4irT. 



g\p') = g 2 (fi) 1 + g 2 m)(3/2-K 2 ) log(fi/ji' 



(2.60) 



All perturbative results on the other hand suffer from the problem of renormalization 
scheme dependence, the more so the stronger the coupling. Having settled for the MS- 
scheme, all of the remaining ambiguity is in the choice of the renormalization scale pf. 
Throughout this paper, we shall choose p = 2nT as our fiducial scale and consider the range 
pf = ttT . . . AnT to test for the scheme dependence of the various approximations.^ 

The leading-order (HTL) result, Eq. (|2.43|) , is simply m/T = m/T = g. For g = g(2nT), 
this is the straight long-dashed line in Fig. [1|. For the different choices p' = TtT and AttT, g 
is instead the function of g(2irT) given by Eq. ( |2.60 ) and m/T is given by the long-dashed 
lines below and above the central one. 

The NLO correction ( p. 45 ) is negative, eventually making the perturbative result for 
m 2 = m 2 + 5m 2 negative, in fact already at moderately large coupling g m 1 (shorter-dashed 
lines in Fig. [I], individually corresponding to p' = ttT,2ttT,4:7tT again). Clearly, using this 
strictly perturbative result would make the thermodynamic potentials fall back to the free 
result at g = tc/3 where m 2 + 5m 2 vanishes, and give rise to tachyonic singularities beyond. 

However, there is no unique "strictly perturbative" result. Defining a NLO mass through 
m = m + 5m would involve 5m = 5m 2 /2m. This would lead to an obvious breakdown of 
perturbation theory only for twice as large values of g, g > 2tt/3 ~ 2, with negative rather 
than imaginary values for m. 



6 In Ref. [22], from which we deviate slightly in taking p = 2ttT rather than p = T as the fiducial 
scale, the scheme dependence of thermal perturbation theory has been studied in the above scalar 
model in great detail with the result that at least at high orders of perturbation theory p! ~ 2irT 
seems to be an optimal renormalization point, corroborating the expectations expressed in Ref. Jl. 
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g(2nT) 

FIG. 2. Comparison of perturbative and HTL-improved approximations to the entropy in the 
large- iV scalar 0(A r )-model. The shaded areas denote the variation under changes of the renor- 
malization scale from p,' = ttT to 47rT. The band marked "HTL" refers to using the leading-order 
(HTL) mass in the 2-loop ^-derivable entropy, "NLA" to using the approximately self-consistent 
NLO mass djgp . Also given are the corresponding results for a naive strictly perturbative NLO 
mass when defined through m 2 or m, respectively. 



But this does not mean that there is no physical content in the NLO effects beyond 
g ~ 1. Rather, the physical content is unnecessarily lost by the restriction to a polynomial 
result for m 2 (or m) which does not preserve the monotonous behavior of m/T as a function 
of g that is observed both at leading order and in the exact result. 

In order to ensure such a monotonous behavior, in Refs. flTBllTBl we have considered the 
simple Pade approximant m 2 + 5m 2 — > m 2 /[l — 5m 2 /m 2 ] = g 2 T 2 /[l + 3g/ir], which already 
achieves a dramatic improvement for g •> 1. An alternative, which is in fact more in the 
spirit of approximate self-consistency, is to return to the approximate gap equation ( |2.42| ) 



m 2 = g 2 T 2 g 2 Tm . 

7T 



(2.61) 



and solve this quadratic equation for m exactly, yielding 

m NLA /T = Jg* + (3g 2 /2n) 2 - 3g 2 /2ir. 



(2.62) 

In what follows, this will be referred to as our "next-to-leading approximation" (NLA) 
for the scalar thermal mass. Also this approximation preserves the propertym of being a 
monotonously growing function of g. For very large g it saturates at winla — y ( 7r /3)T. The 
corresponding results for the various renormalization scale choices are given by the dotted 
lines in Fig. [I], showing a striking improvement over the standard perturbative results also 
for very large coupling. 

With m approximated either by its leading-order (HTL) value m = g(fl)T or by the 
NLA result (|2.61| ), the correspondingly approximated entropy is obtained by evaluating 
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numerically the expression ( |2.38| ). In Fig. this is compared with the strictly perturbative 



expressions for S/S up to and including order g 2 , and g 3 , respectively!] The shaded bands 
indicate the variation of the results with p! = ttT . . . 4nT. Evidently the perturbative 3rd- 
order result fails to be a better approximation than the 2nd-order one for g > 1. The 
semi-perturbatively evaluated HTL result is already an appreciable improvement over the 
2nd-order perturbative result, whereas the NLA follows closely the exact (N — ► oo) result. 
Also shown are the results corresponding to the two "strictly perturbative" NLO mass 
definitions mentioned above when used in the same manner. 



III. QCD: APPROXIMATELY SELF-CONSISTENT RESUMMATIONS 

We turn now to our main case of interest, the QCD plasma. In this section, we shall 
concentrate on a purely gluonic plasma, deferring the addition of quarks to the next section. 
Although the thermodynamic potential in QCD is a gauge independent quantity, in writing 
down its skeleton representation we have to specify a gauge. In formulating the two-loop 
^-derivable approximation we find it convenient to start with the temporal axial gauge. 
While this approximation is by itself gauge dependent, when supplemented by perturbative 
approximations on the generalized gap equation it results in a gauge invariant resummation 
scheme for the entropy. 

A. The skeleton representation of the entropy 

In QCD, the thermodynamic potential is a functional of the full gluon (£)), quark (S), 
and Faddeev-Popov ghost (D gh ) propagators, 



PQ[D,S,D gh ] = -Tr logZT 1 - Tr log^T 1 - Tr logZ^ 1 



1 

2 

-1-TtUD+ TrS5+ Tr U gh D gh + $[£>, S, D gh ], (3.1) 



where Tr now includes traces over color indices, and also over Lorentz and spinor indices 
when applicable. The self-energies for gluons, quarks and ghosts are denoted respectively by 
IT, E and U g h- In Fig. |3], the lowest-order (two-loop) skeleton diagrams for <3> are displayed. 

In gauges which do not break rotational invariance, the gluon propagator at finite tem- 
perature contains up to four different structure functions [ 24| . Only two of them correspond 



to degrees of freedom which are transverse in 4 dimensions; the remaining ones are unphys- 
ical, constrained by a Ward identity ||25|| , and compensated for by the Faddeev-Popov ghost 
degrees of freedom. 

In general, the gluon self-energy 11^ (fc) is a tensor which is not transverse with respect 
to the 4-momentum = (u;,k), but also contains up to 4 structure functions. There are 



7 This plot differs from the corresponding one presented in Ref. IE] in that in the latter the fiducial 



renormalization scale \i = T has been used, so the abscissae are non-linearly related. 
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(b) (c) 




FIG. 3. Diagrams for <3? at 2-loop order in QCD. Wigly, plain, and dotted lines refer respectively 
to gluons, quarks, and ghosts. 

however gauges where ghosts decouple and where as a consequence is strictly transverse^]: 
axial gauges n^A^ = 0, with n M a constant 4-vector. 

A particularly convenient choice appears to be the temporal axial gauge, where coin- 
cides with the rest-frame velocity of the heat bath and thus preserves rotational invariance. 
Ignoring the well-known difficulties with this gauge in the imaginary-time formalism []26|] , the 
temporal axial gauge would lead to great simplifications of the structure of Eq. ( [3.1|) : The 
ghost self-energy U gh vanishes and the ghost propagator does not appear in Secondly, 
there are only two independent structure functions in the gluon self-energy, which can then 
be written as (suppressing the color labels) 

n,>, k ) = L - k -^\ n T (c, k) - k), (3.2) 

n 00 (u;, k) = -U L (u, k), n 0i (u;, k) = — -^U l (uj, k). (3.3) 
With these definitions, the propagator in temporal axial gauge reads 

D™ G (u, k) = L - k M D T (u, k) + ^^ 2 D L (u, k) (3.4) 



where 



-1 -1 

D T[u,k) = — — —, — 7T, D L (u,k) = — r-. (3.5) 

uj — k — \\rp\oj , k) k 2 + Il L {uj,k) 

Note that because Dq^ g = D^q G = 0, only the spatial components Ely of the polarization 
tensor enter Eq. ( |3.1| ) in temporal axial gauge. 

For later use we introduce the following spectral representations: 



8 This property can nevertheless be lost in approximations which do not preserve gauge symmetry; 
cf. the discussion after Eq. ( f3.ll ). 
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D T (u, k) 



00 dA^o Prikp, k) 

-00 2lT kn — 00 



D L (uj,k) = -—+ - . 3.6 

k J -00 2tt k — u) 

Here px and pi are the spectral densities: 

pL.T{ko,k) = Um2lmD L , T (k + ir],k) . (3.7) 

[Note the subtraction performed in the spectral representation of Dl(u, k): this is necessary 
since D L (uj,k) — > — 1/k 2 as |u;| — > 00. At tree-level, = and = 27re(A;o)^(fc 2 ), and 
therefore dJ? } = -1/(uj 2 - k 2 ) and L>[ 0) = -1/k 2 .] 

Concentrating on the gluonic contributions for now and postponing the inclusion of 
fermions to the next section, we obtain in analogy to Eq. (|2.8|) 

Vg/V = N g J-jL* | 2 ( Im log(-o; 2 + k 2 + U T ) - Im U T D T ) 

+ (lmlog(fc 2 + n L )+ lmIl L D L )\ + T$ g [D T ,D L ]/V (3.8) 

where N g is the number of gluons (iV 2 — 1 for SU(iV), i.e. 8 for QCD).[] The entropy of purely 
gluonic QCD can then be written in complete analogy to the derivation of Eq. (|2.12|) as 

S = S T + S L + S' (3.9) 

where 

S T = -2N g I Im log( _^2 + k 2 + Ut) _ Im n T R eJ D T }, (3.10a) 



(2tt) 4 dT 

S L = - N 9 j ^Y^QT i lm ^g{k 2 + Yi L ) + lmU L ReD L }, (3.10b) 



and 

S' 



d{T<&) r d 4 k dn{uj 



dT 



+ N 9 J n ^> (2 Re n r Im D T - Ren L ImD L ). (3.11) 



As in the scalar case, we are interested in the ^-derivable approximation obtained by 
keeping only the two-loop skeletons of Fig. |||. In gauge theories, however, the ^-derivable 
approximations have in general the drawback of violating gauge symmetry, because vertex 
functions are not treated on equal footing with self-energies (in particular, in the two- 
loop approximation to <3> there are no vertex corrections at all). Thus the corresponding 



9 Here we have assumed a principal- value treatment of the factor k 2 jio 2 in Eq. ( |3.4| ) for the contour 
integration. Because this factor is real and positive, it can be dropped from within the imaginary 
part of the logarithm involving 11^. 
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approximation to the polarization tensor IT^ needs not be transverse. Nevertheless, in the 
temporal axial gauge, the previous expressions are not affected by a loss of 4-dimensional 
transversality, because they involve only the spatial components ILj, or equivalently Ut and 
II/. (cf- Eq. (pD). 

Therefore, in this gauge, the property that S' = in the two-loop approximation to 
$ still holds, for the same, essentially combinatorial reasons as in the scalar field theory 
with cubic and quartic interactions of the previous section. In this approximation, the self- 
energies n r , U L and propagators D T , D L are to be determined self- consistently, by solving 
the generalized "gap equations" 

D T l = -u 2 + k 2 + n T [D T , D L ] , D L l = — k 2 — Ii L [D T , D L \ , (3.12) 

i.e., the Dyson equations where U S [D T , D L ] (s = T, L) are the one-loop self-energies built 
out of D T and D L . 

Whereas the entropy expressions ( |3.10| ) themselves are manifestly UV finite, Eqs. Q3.12p 



contain UV divergences which require renormalization. Because of the simple Ward identities 
of axial gauges, (wave function) renormalization of the gluon self-energy at lowest order in 



g contains the correct one- loop coefficient of the beta function [[27],[^| . Beyond lowest order, 
however, it is not clear that the gap equations ( |3.12j ) can be renormalized in a simple manner 
(in contrast to the scalar toy model of Sect. [11 B| ). 

At any rate, in general gauges the 2-loop ^-derivable approximation misses the correct 
perturbative running of the coupling constant. Indeed, the latter is an order-g 4 effect in 
the thermodynamic potentials and is thus beyond the perturbative accuracy of a 2-loop 
^-derivable approximation. 



B. Approximately self-consistent solutions 



Unlike the scalar theory with A0 4 interactions, in QCD the "gap equations" (|3.12|) are 



non-local, which makes their exact solution prohibitively difficult. But in fact, as we have 
just explained, uncertainties concerning gauge symmetry and renormalization beyond order 
g 3 make such a fully self-consistent solution not really desirable. 

For this reason we shall construct approximately self- consistent solutions which main- 
tain equivalence with conventional perturbation theory up to and including order g 3 (the 
maximum perturbative accuracy allowed by two- loop approximations for $), and which 
are manifestly gauge-independent and UV finite. After such approximations — where the 
gluon polarization tensor is transverse and the ghost self-energy (in gauges with ghosts) is 
neglected — , Eqs. (|3.10| ) have the same formal structure in any other gauge, and S' = 
to the same accuracy. We can therefore drop the restriction to the somewhat problematic 
temporal axial gauge. For instance, in the more commonly used Coulomb gauge the gauge 
propagator is given by 

D C G (^, k) = D L (cu, k), D™{u, k) = Uj - D T (u>, k) (3.13) 

and the ghost propagator does not contribute as long as there is no nontrivial ghost self- 
energy; in covariant gauges under the same circumstances, the then propagating ghosts just 
compensate for an additional massless pole that is present in the gluon propagator. 
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With the gauge- independent approximations for Hl,t that we shall obtain from (HTL) 
perturbation theory, the effectively one- loop expressions for the entropy, Eqs. (|3.10| ), con- 
stitute a gauge-invariant approximation to the full entropy. By then computing exactly 
these expressions, we shall obtain a gauge-invariant result which is nonperturbative in the 
coupling g, while being equivalent to conventional resummed perturbation theory up to and 
including order g 3 . 

As generally with thermal field theories ||f|], the perturbative solution of Eqs. ( J3.12 ) 



requires to distinguish between soft (k < gT) and hard (k ~ T) fields, which are dressed 
differently by thermal fluctuations. In (purely gluonic) QCD, and in the Coulomb gauge, 
the hard fields are always transverse, while the soft fields — which may be seen as collective 
excitations of the former |7||| — can be either longitudinal, or transverse. 

Because of the limited phase-space, the leading order (LO) contribution of the soft modes 
to the thermodynamical functions is already of order g 3 [I], so the corresponding self-energies 
are needed only to leading order in g. These are the so-called hard thermal loops 11^ and 
il-r p9|j9|], which in the present formalism appear as the solutions to Eqs. ( |3.12|) to LO in g 
and for soft (k ~ gT) external momenta. They read: 



U L (u, k) = rh 2 D 
1 



U T (u,k) o 



I uj j uj + k 
2k ^ uj — k 
uj 2 ^k 2 



rh 2 D + — ttt" - n 



(3.14) 
(3.15) 



with the Debye mass 



7T 2 jo dk 



The HTL's ( p.!4| ) are manifestly UV finite: they derive from one- loop Feynman graphs, 
but involve only the contribution of the thermal fluctuations in the latter (as opposed to 
the vacuum fluctuations, which are responsible for UV divergences). The corresponding 
propagators are then defined via the Dyson equations ( |3.12| ): 

D t \uj, k) = -uj 2 + k 2 + f[ T (cj, k) , Dl\uj, k) = -k 2 - U l (uj, k) , (3.17) 



Note that, for k ~ gT, the self-energy corrections in Eqs. ( |3.14| )-( |3.17| ) are as important as 



the corresponding tree-level inverse propagators Dq 1 ~ k 2 ~ g 2 T 2 . Thus, at soft momenta, 
the self-energies cannot be expanded out of the HTL-resummed propagators. The HTL 
spectral densities consist of quasiparticle poles at time-like momenta and Landau damping 
cuts for |co>| < k. When k ^> gT, the transverse pole describes the usual single-particle 
excitations (hard transverse gluons), while the additional pole associated to the collective 
longitudinal excitation has exponentially vanishing residue ||31|| . 

For hard, transverse, fields, we need the solution Ii T {k ~ T) of Eqs. ( p,12| ) to leading, 
and next-to-leading order (NLO). This is obtained as: 

Tl T (k ~ T) ~ ILP + 5U T , (3.18) 

where LT^ = n T [D ] ~ g 2 is the bare one-loop self-energy (i.e., the standard one- loop dia- 
grams with tree- level propagators Dq = (Djl , D^) on the internal lines), and 5Ut ~ g 3 T 2 
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FIG. 4. NLO contributions to 5Ht at hard momentum. Thick dashed and wiggly lines with a 
blob represent HTL-resummed longitudinal and transverse propagators, respectively. 

is an effective one-loop self-energy where one of the internal lines is hard (and transverse), 
while the other one is soft (longitudinal or transverse) and dressed by the HTL. Thus, 5Ht 
is the sum of the four diagrams depicted in Fig. f|; these are explicitly computed in App. 

El 

A priori, the one-loop self-energy involves also vacuum fluctuations, and therefore UV 
divergences, which call for renormalization. The UV divergences could be absorbed by 
a wave-function renormalization constant, which drops out from the entropy expressions 
( |3.10| ). As it will turn out presently, only the light-cone limit of 11^ [-Do] will contribute to the 
order of interest. In line with our strategy of restricting to gauge- invariant approximations 
to the self-energy, we shall altogether drop the gauge- dependent vacuum pieces, which in 
fact vanish on the light-cone. 

Because from standard HTL perturbation theory we take UV finite approximations for 
n_L,Tj we shall in fact have no inherent beta function^ prescribing the scale dependence of 
the coupling g. When numerically evaluating the results, we shall simply adopt the standard 
running coupling constant of the MS scheme and consider the resulting renormalization-scale 
dependence of our results as an estimate of our theoretical error (cf. Sect. |II D|) . 



C. Perturbation theory: Lowest orders 

In this and the following subsections, we shall consider the perturbative expansion of 
our master equation for the entropy, Eqs. ( |3.10| ), and recover in the process the standard 
perturbative results up to order g 3 . This is useful not only as a cross check of the various 
approximations, but also as an illustration of the rather non-trivial way that perturbation 
theory gets reorganized by this equation. Moreover, the perturbative expansion will shed 
more light on the physical interpretation of the various terms in Eqs. ( |3 .10 ), and give us 
hints for better approximations to be used in the non-perturbative, numerical calculations 
to come. 



10 A refinement of the present approach which is accurate at and above order g and which has cor- 
rect (lowest-order) coupling constant renormalization would require at least a 3-loop approximation 
to the thermodynamic potentials. 
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The leading-order result is obtained by putting Ut = Hl = in Eqs. Q3.10I) . This is the 
Stefan-Boltzmann entropy of a free gas of massless transverse gluons: 



d 3 k d r . „ t/ ^ n 4tt 2 



-2iV 9 



(2tt) 3 <9T 



Tlog(l - e~ fe/i ) =—N g T A . (3.19) 



Here the retarded prescription (a? — > c<j+ie) is implicit in the first integral, which is evaluated 
with the help of the identities ( gjgg ) and ( ggg ). 



The order-g 2 contribution to the entropy comes also exclusively from hard transverse 
gluons, via one- loop corrections. Specifically, by expanding Eq. ( |3.10a| ) to order g 2 , one 
obtains: 

A7 f d 4 k dn f IiP m 1 

S 2 = -2N a —— — l-lm - , - + ImnV } Re 



9 7 (2vr) 4 <9T | cj 2 - P T u 2 - k 2 

f d A k dn ^ i 2 \ _ 1 



(2tt) 4 dT J cj 2 - A; 2 
= -27riV 3 J-0^^<^ 2 ~ k2 ) Re 4 2) fc) , (3.20) 

where the integral is indeed dominated by hard momenta k ~ T. Note that involves only 
the light-cone projection Relly (u = k) of the one- loop self-energy for (hard) transverse 

(2 s ) 

gluons Ely (uj, k). This projection is a priori UV finite: indeed, gauge symmetry guarantees 

(2) 

that the vacuum contribution to ReU T (u = k) must vanish. Moreover, quite remarkably, 
this projection turns out to be also momentum-independent [[32]. 



U^\u 2 = k 2 ) = g 2 NT 2 /6 = m 2 0O , (3.21) 

and thus defines a (thermal) mass correction, also known as the asymptotic mass. Thus, 
finally, 

5 2 = -N,^f = -^"T, (3.22) 

which is indeed the correct result [|I|. Note also that at leading order the asymptotic mass 
is simply related to the (HTL) Debye mass: m 2 ^ = m 2 D /2. 

It is worth emphasizing that Eq. ( |3.20| ) is the same as the entropy of an ideal gas of 
massive particles (with constant masses equal to m^) when expanded to leading order in 
m 2 ^. As was the case in the scalar model discussed in Sect. II, such a simple identification 
is specific to the entropy, and does not hold for the order-g 2 effect in the pressure. 

In the scalar case we have seen that the HTL-resummed one-loop pressure over-includes 
the LO interaction term by a factor of two. For gluons, Ref. |11[ reported instead a factor 
of three. Inspecting the corresponding calculation reveals that this arises because of an 
incomplete implementation of dimensional regularisation. While in the latter 2 — 2e trans- 
verse polarisations of the gluons are considered, the HTL expressions for have not been 
modified accordingly. However in d ^ 3 spatial dimensions, Eqs. (|3.15|J3.14] ) become 
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d 



' u 2 - e ; 



n, 



m 



D 



1 d fc 2 , 
2 ^2' 1; 2'^' 



(3.23) 



where 



m 



D 



(d-l)g 2 NT 2 (-) d - 3C{d 1)1 
/i 



7r (d+l)/2 



(d>2) 



(3.24) 



as determined by the d-dimensional analog of Eq. Q3.16Q . This gives a real and constant 
11^ = m 2 D = (d — l)^ 2 ^ such that the order-g 2 contribution to the 1-loop HTL-resummed 
pressure Phtl = —\ Tr Xo^Dq 1 + TV) is 



p( 2 ) _ JU 
r HTL — ly 9 



d d+1 k 



(2tt 



I (2+1 



n(uj) Im 



A* 



lu 2 — k 2 



N g m 2 D 



d d+l k 



(2tt 



\d+l 



n(oj) Im 



uj 2 — /c 2 



2 x P 2 (3.25) 



as d — > 3, with dimensional regularization eliminating the quadratic divergence for u — > 
— oo. This is then consistent with momentum cut-off regularization, where d = 3 can be 
kept throughout, after dropping a divergence oc m 2 D A 2 . Presumably, the numerical results 



reported in Ref. [DJ will change significantly when corrected accordingly. 

This sensitivity to (a consistent usage of) regularization schemes is related in fact to 
the UV behavior of HTL-screened perturbation theory; it is not present in our UV-finite 
HTL-resummation of (two-loop) entropy and density. 



D. Perturbation theory: Order g s 

The extraction of the order-g 3 contribution to the entropy in Eq. ( |3.10| ) turns out to be 
more intricate than the standard calculation of the plasmon effect in the pressure M. 



1. The order g 3 in the pressure 

Let us briefly discuss first the plasmon effect in the pressure, as obtained from the skeleton 
representation ( |2.1|) . As explained for the scalar case in Sect. |II Q , the order-g 3 contribution 
to the pressure comes entirely from soft momenta, and reads (cf. Eq. Q2.54 )): 



d 4 k T 

— Im 



(2tt) 4 uj 



log(l + .Dofl) - D fl 



(3.26) 



In QCD, D = (D T ,D L ), tl = (fl T , — fl L ), and a sum over color and polarization states is 
implicit in ( |3.26|) . [Note the minus sign in front of IT^ in these compact notations; this 
reflects our conventions in Eqs. (|3.2|)-(|3~5]).l The integral over uj yields: 



duj 



1TLU 



Im 



log(l + D fl) - flD = log 1 + D Q (uo = 0)fl(u = 0) 



log 1 + 



m 



D 



m 



v 



- f[(u = 0)D (uj = 0) 

(3.27) 
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where the non-vanishing contribution in the second line comes from the longitudinal sector 
alone [33|, since tljj(u} = 0) = m 2 D , while IIr(a; = 0) = 0. Thus, 



Pa 



-N 9 T 



d 3 k 



log 1 + 



m 



D 



k 2 



k 2 



N, 



m 3 D T 



12tt 



(3.28) 



is indeed the stan- 



where the color factor N g = N 2 — 1 has been reintroduced. Eq. 
dard result for P3, generally obtained by summing the ring diagrams in the imaginary-time 
perturbation theory 

The order-g 3 effect in the entropy can be now directly calculated as the total derivative 



of P3 with respect to T. We thus obtain S3 = S3 + S3' , where 



,(b) 



dT 



-N„ 



d A k dn{uj) 



2vr 



dT 



Im 



io g (i + D fi) - nA 
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(3.29a) 



is the derivative at fixed II = fl (recall that the HTL's depend upon the temperature only 
via the Debye mass; cf. Eqs. Q3.14Q and (|3.16|) ), and 



s 



(6) 



dP 3 dm d 



dm D dT 



-N n 



d 4 k 
(2k) 4 



■ n(oj) Im 



dT ( 



Dn 



m 3 D 
Air 



(3.29b) 



This decomposition of S3 is interesting in view of the comparison with the perturbative 
expansion of Eqs. (|3.10|) , to which we now turn. 



2. The order g 3 in the entropy 



chard 
^3 



-N n 



Re5n T (cj = k). 



(3.30) 



Unlike what happens for the pressure, the order-g 3 effects of the hard modes do not 
cancel in Eqs. ( |3.10|) , similarly to what we have observed in the scalar case in Sect. |II C 
Rather, we get a non-zero such contribution by replacing Re — > Re 511^ in Eq. ( |3.20|) 
with 5Ut ~ g 3 T 2 the NLO self-energy correction of hard transverse gluons (cf. Eq. ( |3.18| )) 
This yields: 

d 3 k 1 dn[k) 
(27r) 3 fc dT 

Once again, we need only the light-cone projection of the self-energy of the hard particles. 
What is, however, new as compared to the situation at order g 2 is that Re5IIr(u; = k) is 
not a constant "mass correction", but rather a complicated function of k (see Eqs. (|A16|) 
and ( |A17| )). The calculation of l S 3 iard is deferred to the Appendix, but the final result can 
be anticipated, as we shall see shortly. 

The other contributions of order g 3 come from the soft gluons, which can be longitudinal 

S { l ] + S ( t ] . We have (with n(u) ~ T/u): 



or transverse, and we write iS3° ft 



4 3) 



S 



-N„ 



/^^^{im log(A; 2 + n L ) + lmf[ L ReD L } } 



(3) 
T 



-2N n 



(2nY 
d 4 k 



— { Im 



U! 



log 1 



- k 2 



+ 



n 



T 



k 2 



-Imfl T Re (Dt-DP 



(3.31) 



(3.32) 
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where in the transverse sector, the contribution of order g 2 has been subtracted (cf. 
Eq. ( |3.20| )). More precisely, Eq. (|3.20| ) involves the full one- loop self-energy Ily , while 
the subtracted terms in Eq. (|3.32|) involve only Ur, the HTL. This is nevertheless correct 

(2) 

since IT^ and TIt coincide on the light-cone: 



n T (a/ 



(3.33) 



Ultimately, all the contributions of order g 3 displayed in Eqs. ( |3.3(J| )-( |3.32| ) are soft field 



effects: the quantities S^' and Stf 1 are the LO entropies of the soft gluons, while S^ 3,1 ' 1 is 
the NLO correction to the entropy of the hard gluons induced by their coupling to the soft 
fields (cf. Fig. [|). We expect these three contributions to add to the standard result for the 
plasmon effect in the entropy, namely (cf. Eqs. (ETI 



;(3) 



S 3 = N g m%/ '(3tt). 



(3.34) 



This is verified in the Appendix, where the quantities in Eqs. ( |3.3UD ( p.32| ) are explicitly 
computed, but it can be also understood on the basis of the following argument. 
Eqs. ( |3.30|) -( |3.32|) can be compactly rewritten as 



d A k 

(2%j A to 

d^p dn(p ) 
(2tt) 4 dT 



log(l + D U) - UD - Im n Re (D - D 
Re <ffl Im D , 



(3.35) 



where the sum over color and polarization states is again implicit. The first term within the 
(soft) integral over k is obviously the same as S 3 a \ the temperature derivative of P3 at fixed 
rho (cf. Eq. ( |3.29a| )). It thus remains to show that the other terms in Eq. ( |3.35| ) add to 
Sf\ the piece of the entropy involving the derivative of the Debye mass (cf. Eq. ( |3.29b| )). 
That is, one has to prove the following relation: 



d 4 p dn(p ) 
(2tt) 4 dT 



Re 5U Im Dr 



d A k f dn(oS) 



(2tt) 4 I dT 
+n(u) Im 



ImllRe(Z) - D ) 



dT v 



Dr 



(3.36) 



Eq. ( |3.36| ) is nothing but the general 2-loop identity S' = expanded to the order g 3 . Indeed, 
to order g 3 , Eq. ( |2 . 1 3|) implies: 



<9(T$ 2 



dT 



d 4 k dn(u) 
d J (2vr) 4 dT 



Refllmp - D 



d A p dn(p ) 
(2tt) 4 dT 



Re 511 ImDo, 



(3.37) 



where the first integral is saturated by soft momenta k ~ gT, while the second one is 
dominated by p hard, p ~ T. On the other hand, $3[-D] has the explicit expressionF] 



n This follows by expanding <&[D] in powers of g as follows: $[D] = $>[D ] + (5&[D]/5D)\ Do (D 
D Q )+ ••• = $ 2 + $3+ 
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T$ 3 [D] = ^Tt[U[D ](D-D )] 



which implies: 



d A k 



n(u) Im [ft (D - D c 



(3.38) 



<9(T$ E 



dT 



D 



d 4 k J dn(uj) 
(2^p \~cW~ 



Im [IL(D - D )] + n(uo) Im 



dT 



D-D 



(3.39) 



A comparison of Eqs. (|3.37| ) and ( |3.39| ) immediately leads to Eq. ( |3.36| ). 

Moreover, the soft longitudinal and transverse sectors are decoupled at this order: $3[-D] 
in Eq. (p.38|) is simply the sum of two two-loop diagrams, one with a soft electric gluon, the 
other one with a soft magnetic gluon. The condition S' = can be applied to any of these 
two diagrams separately. It follows that Eq. ( |3.36|) must hold separately in the electric, and 
the magnetic sector. This is explicitly verified in the Appendix, via a lengthy calculation. 
Remarkably, Eq. (|3.36| ) provides a relation between the effects of thermal fluctuations on 
the hard and soft excitations, which are both encoded in the two-loop diagrams for $3: By 
opening up the soft line in $3, one obtains the hard one- loop diagram responsible for the 
HTL II; by opening up one of the hard lines, one gets the effective one-loop diagrams for 
5U displayed in Fig. [|. In the case of the scalar theory, this relation is explicitly verified in 

Eqs. (p37D- (ET59D . 

Let us conclude this subsection on perturbation theory with a comment on the higher- 
order contributions to Sl '■ By inspection of Eq. ( 3.10b ), it is easy to verify that not only 



the LO contribution ~ g discussed above, but also the corrections of order g and g , come 
exclusively from soft momenta. Indeed, one can estimate the contribution of hard momenta 
by expanding the integrand in Eq. ( ft. 10 bp in powers of n^/A; 2 , to obtain: 



Im log(fc 2 + n L ) 



Im n L Re 



imn L lIm(n L ) 
~~ k 2 2 W 



+ 



lmU L ImIT L ReIT L 



A- 2 



k 4 



+ 



k 2 + Tl L 



lmU L Imn L Ren L 



A' 2 



A- 4 



(3.40) 



up to terms of order (Ui/k 2 ) 3 . Remarkably, not only the LO terms, but also the NLO ones, 
of order g A , mutually cancel in the sum of the above equations. Thus, as anticipated, the 
hard modes contribute to Sl only at order g 6 or higher. This shows that our approximation 
scheme is rather insensitive to the unphysical, hard longitudinal modes. This is to be 
contrasted to the direct HTL resummation of the pressure where, to one-loop order, the 
longitudinal sector is sensitive to hard momenta already at order g A , as indicated by the 



presence of UV divergences at this order O . 



3. The HTL entropy 

Since 8Ht(u) = k) is a complicated, non-local function, whose numerical treatment is 
difficult, it is interesting to explore first approximations where 8Ut is set to zero. Specifically, 
let us define the following approximation to the entropy, which is obtained from Eqs. (]3.10|) 
by replacing all propagators and self-energies by their HTL counterparts: 
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Shtl = ~N 9 {2 Im log(-. 2 + k 2 + H T ) 

-2Imfl T ReD T + Im log(A; 2 + f[ L ) + lmfl L ReD L }. (3.41) 

We shall succinctly refer to this as the HTL entropy. Clearly, this is still a non-perturbative 
approximation, since its expansion contains all orders in g. 

A priori, Eq. ( |3.41D is not doing justice to the hard particles, since it uses the HTL 



corrections for both hard and soft momenta (while we know that the HTL's are the LO self- 
energies for soft momenta alone). But it turns out that the order-g 2 effect, which is entirely 
due to the hard fields, is nevertheless correctly reproduced by Eq. Q3.41J) : S^ TL = S 2 . The 



point, as emphasized in Sect. [Ill CJ is that S 2 is sensitive only to the light-cone projection 
of the self-energy, where the HTL is a good LO approximation for the hard modes (cf. 
Eq. 

On the other hand, Shtl contains only a part of the g 3 effect, namely that part which is 
associated with the entropy of soft gluons: indeed, it is obvious that the order-g 3 contribution 
to Eq. ( p.41j ) comes from soft momenta alone, where it coincides with Sf {t = + <Sy , cf . 
Eqs. (|3.31|) - (|3.32|) . Let us therefore study this quantity in more detail (it is the same as the 
first integral in Eq. ( |3.35| )): 

4 3 k = -/^^{lm[log(l + D n)-fl J Do] - Im ft Re (£>-£>„)} 

= Si a) + AS 3 , (3.42) 
where S^ = (dP 3 /dT)\ AD = iS> 3 /4 (cf. Eq. ( |3.29a| )), and the remainder is 



d A k 1 

A4 3) + A5| 3) . (3.43) 



AS3 = Ng J~&!kl {2ImIl T Re( J D T -4 0) ) - Imn L Re (p L - £)f)} 



Remarkably, we have found that the transverse and longitudinal contributions to A1S3 cancel 
within the accuracy that we have reached in a numerical integration of Eq. (|3.43| ) (more than 
8 significant digits). With AS 3 = 0, S^ TL is precisely equal to one fourth of the total g 3 
effect, as it was also the case in the scalar theory with g 2 (j) 4 self-interactions (cf. Sect. |II C| ): 



<j(3) _ <j(3) , <j(3) _ 9P 3 
°HTL = °T + °L — -g^T 



5 3 /4. (3.44a) 



In QCD, however, this property is much more subtle: In the scalar theory, the quantity 
which we call here A1S3 was trivially zero, since Imll = in that case. Here, A1S3 = only 
because a compensation takes place in between the transverse and longitudinal contributions 



12 This is to be contrasted with a direct HTL resummation of the one-loop expression for the 
pressure in QCD along the lines of Ref. [11| — there the HTL corrections contribute throughout the 
hard momentum phase space, while no longer being the right approximation. Instead they give 
rise to artificial UV problems. 
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FIG. 5. The integrand of Eq. ( |3.43| ) after performing the energy integral. The transverse (T) 
and longitudinal (L) contributions do not cancel for each value of k; their sum (full line) vanishes 
only upon integration over all k. 



to Eq. ( p.43| ) , both of which arise from Landau-damping contributions at space-like momenta. 
Moreover, this cancellation occurs only after integrating over all energies and momenta (for 
generic k, the result of the energy integral in Eq. (|3.43|) is non-zero, see Fig. ||). Numerically, 



the contributions to S 



(3) 

HTL 



+ turn out to be 



4 3) + si 3) 



(0.34008738 ... - 0.09008738 . . .)S 3 



(3.44b) 



Let us summarize here the various cancellations which take place at order g 3 in the com- 
plete two-loop entropy: The straightforward perturbative expansion of our master equations 
( pMU| ) leads us to Eqs. ( |3.3U| )-( p.32| ), and thus to the following expression for S3 [recall the 
compact notation introduced after Eq. ( |3.26| )1: 



o osott 1 oh 

O3 — Do + Oq 



c-soft 
^3 



dP, 



S. 



hard 



dT 
dP 3 



drriD 



dm D dT 
dP 3 di%D 
dm D dT 



- N„ 



d^k dn(u>) 



(2 



7T 



dT 



ImftRe (D - D c 



- Ad? } - A6i 3) 



(3.45) 
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In these equations, AS^ and AS^ have been defined in Eq. ( |3.43|) , and the second line 
in the above expression for S^' 1 follows either by using S' = (cf. Sect. III.D.2) 
by explicitly computing Eq. (|3.30| ) within HTL-resummed perturbation theory (cf. 
Appendix; see especially Eqs. (|A22|) and (|A23|) there). Furthermore, by construction, S 



:(3) 



, or 
the 

(3) 

HTL 



is the same as Sf ft . 



According to these equations, the quantities AS L and A<S 



:(3) 



cancel in Sf 



independently in the longitudinal and transverse sectors, thus yielding the correct result 
for 1S3, cf. Eqs. ( |3.29a|) and ( |3.29b|) . This is what we have been able to prove analytically 
(cf. Sect. III.D.2 and the Appendix). On the other hand, we have found numerically that 
A<S 3 = AS L 3 ^ + AS^ = 0, so that the actual results for <S§ oft and S^ avd are even simpler: 



c-soft 
^3 



<9A 



dT 



chard 
^3 



dP 3 dm D 



drhn dT 



(3.46) 



At this stage, we have no fundamental understanding of the "sum rule" A1S3 = 0. But this 
serendipitous result will have important consequences in practice, as we shall see below, 
because it determines the magnitude of ,5>3 ard to be 3/4 of S 3 , as was the case in the simple 
scalar model of Sect. [II C , while being an incomparably more complicated expression than 

A full numerical evaluation of the HTL entropy, non-perturbative in g, will be presented 
in Sect. |V A] below, and estimates of the effects of including 6~Ut ~ g 3 T 2 in Sect. |V B|. 



IV. QCD: ADDING THE FERMIONS 

It is now straightforward to add fermions to our theory. We consider Nf flavors of 
massless fermions with equal chemical potential fi; we choose fi > 0, which corresponds 
to an excess of fermions over antifermions for all flavors. Adding the fermions will have 
two effects: first, this will modify the parameters of the gluonic sector, namely the Debye 
mass m 2 D , and therefore also the asymptotic mass m 2 ^ = rh 2 D /2; second, there will be new 
contributions to the entropy. In addition, at finite fi, there is a new thermodynamic function 
of interest, namely the density Af, which shares many of the interesting properties found for 
S. 

The full (leading-order) Debye mass in the QGP reads : 



m 



D 




= (2N + N,)^-+N f ^. (4.1) 

We have introduced here the statistical distribution functions for fermions (/+) and an- 
tifermions (/-), 

and we have used the following integral: 

Jdkk(f + (k) + f4k)) = ^f + Y- ( 4 - 3 ) 



33 



A. Entropy and density from the skeleton expansion 



To construct the fermion contribution to the entropy, let us return to the full skeleton 
representation of the thermodynamic potential (in a ghost-free gauge) and add fermions to 
it. This becomes 

(3Q[D, S] = ^Tr log IT 1 - ^TrUD- Tr log S' 1 + Tr £S + $[£>, S], (4.4) 

where S and S denote respectively the fermion propagator and self-energy, and the sum over 
the gluon polarization states (two transverse and one longitudinal) is implicit. $[D,S] is 
the sum of the 2-particle-irreducible "skeleton" diagrams constructed out of the propagators 
D and S. Below, we shall be mainly interested in the 2-loop approximation to <$>[D,S], 
where the only new diagram is the one represented in Fig. ||]d. The self-energies £ and II in 
Eq. (O) are themselves functionals of the propagators, defined as 

The self-consistent propagators D and S are obtained by solving the Dyson equations 

zr 1 = Dq 1 + n, s- 1 = + s. (4.6) 

Then, the functional $[-D, S] is stationary under variations of D and S around the solutions 
to Eqs. M : 

5Q[D,S]/5S = 0, 5Q[D, S]/5D = 0. (4.7) 

The entropy S(T, /x) and the density Af(T, fi) are obtained as the derivatives of the thermo- 
dynamic potential with respect to the temperature, and the chemical potential, respectively: 



S 



d(Q/V) Kf d(Q/V) 



dT 



<9/i 



(4.8) 



Because of the stationarity property (|4.7|), we can ignore the T and /i dependences of the 
spectral densities of the propagators when differentiating <&[D,S]. That is, we have to 
differentiate only the statistical factors n(to) = l/(e |8a; — 1) and f(to) = l/^e^ - ^ + 1) which 
arise after performing the Matsubara sums in Eq. ( |4.4|) . This yields, for the entropy, 

= _ disi/v) 



dT 



^St + Sj + S*, (4.9) 



where Sb = St + <Sl is the purely gluonic part of the entropy, as shown in Eqs. ( |3.9| )- (|3.10b| ), 
Sf is the corresponding fermionic piece, which reads (the trace below refers to Dirac indices) 

Sf = -2 j** tr { Im M70S- 1 ) - Im ( 7o S) Re (S 7 o)} , (4.10) 

and 
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S' EE 



d(T$) 



dT 



D,S 



d 4 k 



dn(uj) „ df(uj) 

' RenimD + 2-^^tr 



<9T 



<9T 



Re( 7 oS)Im(5 7 o) (4.11) 



has the important property to vanish at 2- loop order That is, 5 ~ + iSy to the order 
of interest. 

The corresponding expression for the density is obtained by replacing (df /dT) — > 
(df /dfi) in all the formulae above. This gives M = TV/ + A/"', with M' = in the 2-loop 
approximation. Thus, to the order of interest, 



d 4 k df(u) 
(2vr) 4 dn 



tr { Im \og{ l0 S- 1 ) - Im ( To E) Re (S 7o )} 



(4.12) 



For simplicity, all the previous formulae have been written for only one fermionic degree 
of freedom; the corresponding formulae for iV colors and Nf flavors can be obtained by 
multiplying the fermionic contributions above by NNf. 

Note finally the following Maxwell relations, 



dS 

dfj, 



dAf 



dT 



(4.13) 



which express the equality of the mixed, second order derivatives of the thermodynamic 
potential. In our subsequent, self-consistent construction of S and M, these relations will 
be satisfied at the same order as the requirement of self-consistency. 



B. The structure of the fermion propagator 

In the previous formulae we have always associated a factor of 70 with the fermion 
propagator and self-energy. This was possible since 7q = 1 and det7o = 1; it is also 
convenient since, e.g., S> = 70S70, and it is preferable to work with hermitian Dirac matrices. 

In order to compute the Dirac traces in Eqs. (^4.10|) -( f4.12|) , it is useful to recall the 
structure of the fermion propagator at finite temperature and density: The most general 
form of the self-energy E which is compatible with the rotational and chiral symmetries is: 

k) = a(w, k) 7 + b(u, k)k ■ 7. (4.14) 

(For a massive fermion, this would also include a mass correction, i.e., E = a{uo 1 k) r y° + 
b(u>, k)k ■ 7 + c{uj, k).) This can be rewritten as: 

7oE(cu,k) = E + (cj,A;)A + (k) - E_(c^, fc)A_(k), (4.15) 

where E-j-(a?, k) = b(u, k) ± a(u, k), and the spin matrices 

A ± (k)EE A + + A_ = l, 

A^ = A ± , A + A_ = A_A + = 0, trA ± = 2, (4.16) 

project onto spinors whose chirality is equal (A+), or opposite (A_), to their helicity. Dyson's 
equation S^ 1 — — }t + E then implies: 
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-f S-\u;,k) = A-\u,k)A + + Az\uj,k)A-, (4.17) 
with A^ 1 = -[u T(k + E±)]. This is trivially inverted to yield the fermion propagator: 

S-fo(co,k) = A + (cu,k)A + + A_(u,k)A_. (4.18) 



The presence of the projection operators A± in Eqs. (|4.15| ), (|4.17| ) and ( [4.18|) allows one 



to easily compute the Dirac traces in Eqs. ( [4.1 0| ) and ( |4.12j ), and thus obtain: 



= "4 {Im log A; 1 + Im ^(-Al 1 ) + 

- Im£ + ReA + + ImXLReA-}. (4.19) 

The corresponding expression for Af is obtained by replacing (df/dT) — > (df/dfi) in the 
equation above. 

C. Perturbation theory for Sf : order g 2 



Eq. ( }4.19|) will be now supplied with certain approximations for the quark self-energies 
E±. As before, we aim at reproducing the results of perturbation theory up to order g 3 . 
This will be achieved by approximations analogous to those employed for the gluons, namely 
the HTL approximation E±, supplemented by the NLO correction 5'E± to the hard fermion 
self-energy on the light cone. 

Note, however, an important difference with respect to the gluon case: unlike the soft 
gluons, which contribute to the entropy already at order g 3 , the soft fermions contribute 
only at order g A or higher, because their contribution is not enhanced by the statistics. Nev- 
ertheless, in our numerical calculation below, we shall carefully include the contribution of 
the soft fermions, appropriately dressed by the HTL. This is in line with our general strategy 
of constructing non-perturbative approximations for the entropy (or other thermodynamic 
quantities) which include as much as possible the dominant collective effects in the plasma. 

In the HTL approximation, the fermion self-energies read as follows [|30]|§ : 




W )^ |1 -^ 1 „ g ^|), (4,0) 

where M 2 is the plasma frequency for fermions, i.e., the frequency of long- wavelength {k — > 0) 
fermionic excitations (C/ = (A^ 2 — 1)/2N): 

= g -^f^k(2n(k) + f+ (k) + Uk)) = ^(T 2 + g). (4.21) 

We are now in position to evaluate the fermionic entropy and density up to order g 2 : 
To zeroth order, i.e., for an ideal gas of massless fermions at temperature T and chemi- 
cal potential fi, we obtain the well known results fj] (the color-flavor factor NNf is here 
reintroduced): 
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c(0) 
6 / 



NNt 



' 7tt 2 T 3 /x 2 T n 
45 + 3 . 



AA (0) = NN-, 



T 2 + 



71- 



(4.22) 



The correction of order g 2 involves the fermion self-energies to one loop order, E±^ : 

Sf/NN f = -4 % { Re E? Im -=L - Re eL 2) Im -=L} 

df(-k) 



(2vr)< 
" (2vr) 3 \ <9T 



ReE^V = k) 



dT 



ReEL 2) (a; 



-k) 



(4.23) 



As in the gluon case (cf. Eq. ( 3.20Q ), the correction of order g 2 is sensitive only to the 
light-cone projection of the one-loop self-energy, which is correctly reproduced by the HTL 
approximation ( 4.20| ) [32| . That is, 



ReEi 2) (cj = ±k) 



M 2 

Z ± (LU = ±k) = — . 



(4.24) 



Eqs. (|4.18|) and ( f4.24|) show that, to order g 2 , the hard fermions (or antifermions) propagate 
as massive particles, with dispersion relation ej, = k 2 + 2M 2 . This identifies the fermionic 
asymptotic mass as = 2M 2 . By also using the properties f(k) = f+(k) and f{—k) = 
1 - /_(&) (cf. Eq. (£l|)), together with Eq. fl4~3| ), we finally deduce 



<5<®WN - M 2 d 
S f / NN f ~ ~-^fff 



7f 2 T 2 



IS 
2 



M 2 T 



6 



(4.25) 



The leading-order correction to the density Afi 2 ^ is obtained similarly: 



n 2 T 2 



EL 
2 



[iM 2 



2 

oo 



7T 



2tt 2 



(4.26) 



The above results for and Aj 2 '*, together with the previous ones for scalars, Eq. Q2.57| ), 
or gluons, Eq. ( |3.22|) , can be generalized to the following, remarkably simple, formulae, which 
hold for an arbitrary field theory involving massless bosons (with zero chemical potentials) 
and fermions: 



S- 2 



-t E 

I B 



m 



ooB 



12 



+ E 

F 



24 



AA 



57T 



ooF- 



(4.27) 



Here the sums run over all the bosonic (B) and fermionic (F) degrees of freedom (e.g. 4 
for each Dirac fermion), which are allowed to have different asymptotic masses and, in the 
case of fermions, different chemical potentials. According to Eq. ( |4.27| ), the leading-order 
interaction term in the entropy as well as in the density has a very simple physical origin: 
it is entirely due to the thermal masses acquired by the hard plasma particles, i.e., directly 
given by the spectral properties of the dominant degrees of freedom. 

To conclude this discussion of the order g 2 , let us summarize here the respective contri- 
butions to entropy (S2 = + <Sf^) and density (A/2) in hot SU(iV) gauge theory with Nf 
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quark flavors: these follow from Eqs. ( |3.22|) , ( |4.25| ), ( |4.26|) (with the thermal masses ( |4~T 
and (|4.21| )), and read: 



g 2 N g T (AN + 5N f _ 2 , 3N f ,1 g 2 »N g N f ( ' 2 , /i 



P2 -"^2"\ 18 T + ^ /iT + 2^ /i J- (428) 

In writing these equations, we have also added the corresponding expression of the pressure 
(P2), as taken from Ref. JTJ. Clearly, our above results for S2 and A/2 are consistent with 
this expression for P 2 : S2 = dPijcfT, A/2 = dPijd\i. 



D. Perturbation theory for 5/ : order g 3 



Unlike the g corrections in Eq. ( |4.28| ), — which apply to the whole area of the \x— T 
plane where the coupling constant is small (i.e., such that max(/i, T) is much larger than 
^■qcd) — j the corrections of order g 3 that we shall discuss now apply only to the high 
temperature regime^] T 3> rhp. This restriction is obvious in the imaginary time formulation 
of thermal perturbation theory, where the effects of order g 3 arise entirely from the sector 
with zero Matsubara frequency H . In the present calculation, these effects are obtained by 
approximating n(k) ~ T/k for k ~ rho-, which is valid provided rfi£, T. Assuming this 
condition to be satisfied, we shall now show how the "plasmon effect" arises in our formalism 
when the fermions are also included. This is similar to the previous discussion of the pure 
glue case (cf. Sect. [Ill D| ), so we shall indicate here only the relevant differences. 

There are two types of contributions of order g 3 to the entropy: (i) the direct contribution 
of the soft gluons, «S| oft = +5^ 3) , which is still given by Eqs. (|3~31p and ( ^32|) , and (it) the 



NLO correction 53 ard to the entropy of the hard particles, which now includes contributions 
from both transverse gluons and fermions, via the NLO corrections to the corresponding 
self-energies on the light cone (cf. Eq. ( |3.30| ) and ([4.23|) ): 



(2nf 

+ 2NN f (^RRe5Z + (uj = k) + Re 5S_(cu = -*)) }. (4.29) 

The diagrams pertinent to 5I1t have been shown in Fig. |4j. The corresponding diagrams for 
5H± are similar, and are displayed in Fig. || Their evaluation proceeds along the same lines, 
and is briefly discussed in App. 0. Let us summarize here the final results: 

As in the pure glue case, it can be verified that there is no net contribution from the 
soft transverse gluons: the direct contribution in Eq. (|3.30|) is precisely cancelled by 
the corresponding contributions to the self-energies of the hard particles, 511^ and <5£j_ (cf. 



13 If fi = 0, then mc~ gT, and this condition is equivalent to weak coupling; for \i > 0, however, 
there is a new scale in the problem, and the high-T condition becomes an independent condition. 
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FIG. 6. NLO contributions to 5T, at hard momentum. Thick dashed and wiggly lines with a 
blob represent HTL-resummed longitudinal and transverse propagators, respectively. 

Figs. H and ||). As expected, the whole contribution of order g 3 comes from soft longitudinal 
gluons (either directly, via «S| , or indirectly, via their contribution to iS3 ard ), and reads: 

N g m 3 D dm 2 D N g m D N g 
12n dT 8vr ~ 12n 

where we have introduced the notation 



S 3 = + T ^ = (rkl + 3m D m»), (4.30) 

tl 

= m 2 T + mj, m| = (2N + N f )^- } m 2 =N f ?£-, (4.31) 



so that TdTm 2 D = 2m\. Note that, formally, Eq. ( |4.3U| ) would predict a no n- vanishing en- 
tropy in the zero temperature limit, coming from the term m\\ this is, however, wrong, since, 
as already mentioned, this expression has been obtained on the basis of a high temperature 
expansion and cannot be extrapolated to small temperatures. 

Still as in the pure glue case, the two terms in the r.h.s. of Eq. (|4.30f) are the same 
as <S>!° ft and i5>3 ard , respectively, because of the "sum rule" AS 3 = 0. (Cf. the discussion 
in Sect. [Ill D 3| ; the arguments leading to Eq. ( |3.46|) are not changed by the addition of 



fermions, since they hold for any value m D of the Debye mass.) The only difference with 
respect Sect. |III D 3| is that, for fi ^ 0, the two terms in Eq. ( |4.30| ) are no longer equal to 
1/4 and, respectively, 3/4 of the total result (compare to Eq. ( |3.44j )); indeed, the identity 
TdTth 2 D = 2m 2 D is valid only at /x = 0. 

Consider now the order-g 3 effect in the quark density: since soft fermions do not con- 
tribute to order g 3 , the only such contribution comes from the NLO corrections 5H± to the 
hard fermion self-energies. This is calculated explicitly in App. |A] along the same lines as 
for the entropy (cf. Eqs. ( |A2S| ) and ( |A29| )) with the result 

Mz = N ' T r_ Dml - ^W. (4.32, 



47T/1 8tt 

The previous expressions for A/3 and 1S3 verify the Maxwell relation 

dS 3 _ cW 3 _ g 2 N g N f fi(m 2 D + m 2 T ) 



(4.33) 



d/j, dT 8tt 3 rho 

which is as expected, since our calculational scheme has preserved self-consistency up to 
order g 3 . These are also consistent with the well-known result for the sum of the ring 
diagrams P 3 = N g Tm 3 D /\2ir . As emphasized already, this result is valid only for high 
enough temperatures, T>mo. In the opposite limit T = 0, it is well known @,^] that 
the sum of the ring diagrams gives a result P r i ng ~ g 4 {i 4 \ogg. 
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V. QCD: NUMERICAL EVALUATIONS 



In the following, we shall turn to a full numerical evaluation of the entropy and the 
density in the approximation S' = = J\f' when further approximated, firstly by the HTL 
approximation (cf. Sect. [Ill D 3|) , secondly by also including NLO corrections to the self- 
energy of hard excitations. 



A. HTL/HDL approximation 

We have seen that the HTL approximation (or in the case of T = and high \x the 
hard-dense-loop [HDL] approximation) is sufficient for a correct leading-order interaction 
term in entropy and/or density — in contrast to a direct HTL approximation of the one- loop 
pressure. On the other hand, the so-called plasmon effect of order g 3 is included only partly, 
namely only in the form of "direct" contributions from soft modes; a (larger) "indirect" 
contribution is due to NLO corrections to the self-energy of hard particles on the light-cone 
as given by standard HTL perturbation theory. 

Since we have found in our scalar toy model of Sect. |11 C| that already the HTL approx- 
imation in the entropy expression with S' = is an improvement over the leading-order 
perturbative result, we shall first concentrate on numerically including all the higher-order 
effects of HTL/HDL propagators in entropy and density. 

Concerning the contributions of the gluonic quasiparticles, the task is to evaluate 
Eq. ( |3.41|) without expanding out the integrand in powers of rhu/T oc g. 

SHThiT^rriE,) involves two physically distinct contributions. One corresponds to the 
transverse and longitudinal gluonic quasiparticle poles, 



cQP __/v 



^[2Tlog(l - e-^) +Tl0g \ e _ e _ t/T ], (5.1) 



where only the explicit T dependences are to be differentiated, and not those implicit in the 
HTL dispersion laws u>T{k) and u^k). The latter are given by the solutions of lo\ — k 2 = 
n T (cu T , k) and k 2 = -fl L (uj L , k) with fl L and fl T given by Eqs. (|3.14|J3.15| ). 



Secondly, there are contributions associated with the continuum part of the spectral 
weights. These read 

30 k 2 dk f dn{ijj) 





-2Imfl T Re [k 2 - to 2 + fir]" 1 + arg[k 2 + fl L ] - Imfl L Re [k 2 + UlV 1 }- (5.2) 



Both the Stefan-Boltzmann part S$b and the standard perturbative g 2 -contribution S2 
of Eq. (|3.22|) are contained in the first term of Eq. ( |5.1| ); all the other terms in Eqs. (|5.1J ),( [5.2|) 
are of order g 3 in a small-g expansion. However, if such an expansion were truncated beyond 
order g 3 , the resulting entropy would be a function of g that initially decreases with g, but 
eventually grows without bound to values larger than S$b (dashed line in Fig. [7]). 

On the other hand, the full numerical result for the HTL entropy (full line in Fig. |7|) 
turns out to be a monotonously decreasing function of rhj^/T. In the case of Eq. (|5.1| ), the 
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FIG. 7. The HTL entropy per gluonic degree of freedom normalized to its Stefan-Boltzmann 
value as a function of the Debye mass rhu(T, fj,)/T. The full line gives the complete numeri- 



cal result corresponding to Eq. fl3,41| ); the dashed line corresponds to the perturbative result to 
order (m,£>/T) 3 ~ g 3 . The dotted line gives the entropy for scalar degrees of freedom with mo- 
mentum-independent mass to = TOqc 
dash-dotted line. 



m.£)/\/2; its perturbative approximant is given by the 



numerical evaluation involves solving first numerically the transcendental equations for ujx{k) 
and UL(k), and a numerical integration, in which it is advisable to separate off the Stefan- 
Boltzmann value through the replacement log[l — e~^ T<yk ^ T ] — > log[(l — e~^ T<yk ^ T /{I — e~ fe / T )]; 
Eq. ( |5.2p requires two successive numerical integrations. 

It is interesting to compare the rather complicated expression Shtl with the simple 
scalar expression Eq. Q2.38|) of the entropy of an ideal gas of massive bosons, 2N g S Q (m), 
which is basically what is considered in the simple massive quasiparticle models of Ref. [|5],^| . 
If in the latter the boson masses are identified with the asymptotic mass of the gluons, 
m = moo = m£)/v^2, then this reproduces the correct leading-order interaction term in the 
entropy. The plasmon effect (i.e. the order-g 3 contribution) is included only partially, but 
not as 1/4 of the complete plasmon effect, but as 1/(4-^/2)- This is because a constant 
thermal mass equal to its asymptotic value underestimates the Debye mass by a factor of 
1/V2 and therefore the plasmon effect by (l/v^) 3 , which is only partially compensated by 
now having 2N g degrees of freedom exhibiting the analog of Debye screening instead of only 
the N g longitudinal ones. 

Numerically, 2N g S (m QO ) reproduces the HTL entropy very accurately (within < 0.1%) 
up to mo ~ T. For larger values of Ad, the HTL entropy leads to significantly larger 
deviations from S$b- This latter fact is somewhat surprising since the plasmon effect, which 
in the HTL entropy is 30% greater than in the simple massive quasiparticle entropy, always 
counteracts the leading-order interaction contribution, as can be seen from the perturbative 
approximant of 2NgSo(m 00 ) (dash-dotted line in Fig. |7|) and that of the HTL entropy (dashed 
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[S HTL - 2N g S (m oo )]/\2N g S {m oo ) - S SB 

0.1 

0.05 


-0.05 
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m D /T 

FIG. 8. Relative deviation of the HTL entropy from that of a gas of massive bosons with 
(constant) mass (full line). The relative deviation of just the transverse quasiparticle contri- 
bution is given by the uppermost dash-dotted line; the transverse Landau-damping contribution 
is given by the lower dash-dotted line. The short-dashed line gives the longitudinal quasiparticle 
contribution; the long-dashed line the longitudinal Landau-damping one. 

line) through order g 3 . This is partly due to the fact that Shtl contains also a term 
~ g 4 \og(c/g), which is not present in the simple massive quasiparticle entropy 2N g S (m oo ) . 

Inspecting in more detail the numerical deviation of the HTL entropy from that of a 
massive gas of bosons, one finds that the quasiparticle contribution from the transverse 
modes, which is always the dominant contribution to the entropy, by itself is always above 
2N g S (m oo ). The transverse Landau-damping contribution is also positive, but relatively 
smaller. On the other hand, both the longitudinal quasiparticle and Landau-damping con- 
tributions are negative, resulting in a small net deviation from the simple massive boson 
entropy for small values of rho/T. When normalized to the deviation of 2N g So(m 00 ) from 
the Stefan-Boltzmann result, the deviation of Shtl from 2N g S (m oo ) is less than about 
+ 1% for rh D /T < 0.739, while negative and rapidly growing for larger values of m D /T, as 
shown in Fig. ||. 

The formulae for the fermionic contributions to the entropy are quite analogous to the 
gluonic contributions. They read Sj^htl = SfJj TL + Sj' 1 ] ITL with 

Sf HTL = NN f /^^{T log(l + c -^«-^) 
o 

1 _|_ e -[w-(ft)-/i]/T 

+Tlog i + e -(fc-M)/r + 0* ( 5 - 3 ) 

where again only the explicit T dependences are to be differentiated and not those implicit 
in the dispersion laws uj + (k) and oo^(k) of the fermionic quasiparticles, which are given by 
the solutions of u± — ±[p + S±(cu±, k)\ with E ± given by Eq. (|4.20|) . 
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FIG. 9. The HTL entropy per quark degree of freedom at \jl = normalized to its free value 
as a function of the fermionic plasma frequency M/T. The full line gives the complete numerical 
result corresponding to Eq. ( |4,10| ) in the HTL approximation; the dashed line corresponds to the 
perturbative result to order (M/T) 2 ~ g 2 ; the dotted line gives the entropy for a fermionic degree 
of freedom with momentum-independent mass M = M^ = y/2M, which has the same perturbative 
approximant to order g 2 . 



The fermionic Landau- damping contribution to the entropy is 



°f,HTL 



-NNf 



k 2 dk 



duo 



df+{u>) , 



|arg[A; - uj + £_ 



dT dT 

— Im £+(tu, k) Re [k — uo + 
+ arg[k + uo + t-(u,k)] - Im £-(u, k) Re [k + u + t-(u, k)}' 1 } 



(5.4) 



In the case of the gluonic contributions to the HTL entropy, there was no difference 
between vanishing and non-zero chemical potential other than the resulting different value 
of ifiD-, which depends on \x according to Eq. ( fOj) . For the quark contributions to the 
entropy, the chemical potential enters both explicitly through the Fermi-Dirac distribution 
function / and through the magnitude of the fermionic plasma frequency M. 

In Fig. |S] the results of a numerical evaluation of the fermionic contribution to the HTL 
entropy normalized to its free value is given as a function of M/T for /j, = 0. When 
compared with the free entropy of simple massive fermions of mass M = Moo = V2M, one 
finds that the HTL entropy exceeds the latter by at most +1.2% for M/T pa 1, coincides 
with it at M/T pa 1.39, and becomes lower for larger M/T.0 On the other hand, the 
strictly perturbative result up to order g 2 is significantly lower, but compared to the gluonic 



14 A 

gain, this good agreement requires all quasiparticle and Landau-damping contributions to- 
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contribution the discrepancy is much smaller because there is no (direct) plasmon effect in 
the fermionic contributions [all order g 3 contributions eventually arise from NLO corrections 

to Moo]. 

Turning now to the quark density, its quasiparticle and Landau-damping contributions 
are obtained by replacing d/dT in the above formulae ( |5.3| , |5.4| ) by d/dfi. 

In the limit T — > 0, the resulting expressions can be simplified to read (for u, > 0) 



HDL 



T=0 



NNf 



k 2 dk 

— —\9(li - uJ+(k)) - 9(u-(k) - /i)]. 

7T Z 



(5.5) 



and 



HDL 



T=0 



— 3-{arg[A;-^ + E + (^A;)] 



— ImS + (//, k) Re[k — /i + £+(/•*> k)] 1 
+ aig[k + fi + 'E-(ji,k)] - ImS_(/i, k)Re[k + u, + £_(//, fc)] -1 }- 

For /x > M, the quasiparticle contribution (|5.5|) can be more explicitly written as 



(5.6) 



1 



3tt 2 



m 3 -^(m) 



(5.7) 



— - u 3 - fci(w) 

t=o,h>m 3vr 2 3vr 2 + vp; J 
where k±(u) is the solution of u)±{k±) = //. 

The first term on the right-hand side of Eq. (|5.7Q represents the free contribution of 
one massless Dirac fermion, the two bracketed terms are the corrections from the nontrivial 
dispersion laws of the two fermionic quasiparticle branches.^ 

For comparison, the fermion density of a free massive Dirac fermion with mass M is 
given by 

_ f _i^( At 2_ M 2> ) 3/2 for fi> M 

t=o ~ 1 for jj, < M 



(5.8) 



Identifying^] M = = \[2M gives the correct leading-order interaction term of order g 2 , 
while leading to somewhat larger values for TV than the perturbative order-g 2 result for all 
M/ p. 



gether; for instance, the normal (+) quasiparticle pole contributions alone would give deviations 
which go up to about +7% for the range of M/T considered. 



15 Because of the "plasmino dip", Eq. ( |5.7| ) becomes more complicated for fi < M, but this case 
corresponds to much too strong coupling to be taken seriously anyway. 



16 Occasionally [12], in simple quasiparticle models of the pressure of fermions at high density the 
identification M = M is made. This happens to give the correct leading-order interaction term of 
order M 2 /^ 2 ~ g 2 there, but only because of compensating errors. At high densities the mass of 
quasiparticles at the Fermi surface is actually = 2M 2 , but in the pressure the leading-order 
interaction term is over-included by precisely a factor 2 when considering only the expression for 
free particles and inserting a constant thermal mass. 
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FIG. 10. The HDL quark density per quark degree of freedom at T = normalized to its 
free value as a function of the fermionic plasma frequency M//i. The full line gives the complete 
numerical result corresponding to Eq. (|4.12|) in the HDL approximation; the dashed line corresponds 
to the perturbative result to order (M//i) 2 ~ g 2 ; the dotted line gives the density for a fermionic 
degree of freedom with momentum-independent mass M = = y2M, which has the same 

perturbative approximant to order g 2 , and vanishes for M > /i, i.e. M/fi > l/\/2. 



In Fig. [TO] the numerical result for Mhdl at T = is plotted for M/fi up to l/y/2, where 
the fermion density of Eq. ( |5.8| ), displayed by the dotted line, vanishes. The HDL result, 
which is given by the full line, is seen to drop to zero almost at the same ratio, to wit, M / fi m 
0.69264. Beyond this point the result becomes negative, showing that the approximation is 
breaking down at such high values of M / \i. [Note that, since M 2 = g 2 [x 2 /&ir 2 for N = 3 and 
T = (cf. Eq. ( [4.21] )), M / ji « 0.69 corresponds to a relatively large coupling g « 5.3.] 

For comparison, the strictly perturbative result to order g 2 is given by the dashed line 
in Fig. [To], which is seen to approach zero faster than the HDL density as well as that of a 
simple massive quasiparticle. 



B. Estimate of NLO contributions 

As we have discussed at length in the previous sections, the HTL approximation in the 
entropy contains only part of the plasmon effect, a different source of order g 3 contributions 
comes from NLO corrections to the gluonic and fermionic self-energies at hard momenta 
on the light-cone as given by Eq. ( |4.29|) . From the result (|3.44|) we know that this NLO 



contribution corresponds precisely to the second term of the right-hand side of Eq. ( |4.30|) . 

In the case of the density, it is clear from the absence of a bosonic distribution function 
in Eq. Q4.12 ) that M in the HTL/HDL approximation does not contain any g 3 contribution, 



so all of A/3 as given by Eq. Q4.32j ) arises from the NLO correction to the quark self-energy 
at hard momenta on the light-cone. 
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As can be seen e.g. from Eqs. ( |A16|) and (|A17 ), the NLO self-energy corrections are 
complicated and nonlocal quantities. Even when evaluated on the light-cone, they do not 
simply give a constant correction to the asymptotic mass, but a nontrivial function of the 
(hard) momentum. In fact, there are even contributions of the form g 2 rhr)P, which grow 
larger than g 3 T 2 for p ^> T, eventually causing a break-down of standard HTL perturbation 
theory, but fortunately such contributions are irrelevant thanks to the fact that n(p) shuts 
off exponentially then. 

Because a full inclusion of the NLO self-energy corrections is rather complicated and 
computationally expensive, and because in the applications below the magnitude of the 
NLO corrections, when treated along the lines of the scalar toy model in Sect. IIP , turns 
out to be comparatively small, we shall in the following consider the approximation of an 
effective constant NLO asymptotic mass. The complete evaluation of 5U and ££, which 
involves a number of technical intricacies, will be reserved for a separate publication. Their 
eventual numerical effects on the thermodynamic potentials is work in progress, though we 
do not expect them to deviate too much from the estimates derived in this subsection. 

From the requirement that a replacement of m 2 ^ and in Eqs. ( |3.22| ) , ( |4.25| ) , and ( f4.26| ) 
by effective constant (i.e. averaged) corrections 8m\ 
(cf. Eqs. ( [4 . 3 0| ) and ( |4.32| ), respectively), we have 



and 5 equals S 3 hard and W 3 hard = A/" 3 



^mlT - l -NN f 5MlT 



2vr 2 



NNfdMlv 



—N g m D rriT 
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(5.9) 
(5.10) 



with rhjj, 
solution 



m T , 



and m M as defined in Eq. ( |4.31| ). This has the remarkably simple unique 



5m 



-— g 2 NTm D , 5Ml 

ZlT 



Q CfTrhn, 

2tt 1 



(5.11) 



where in the latter Cf = N g /(2N). Indeed, in Eq. ( |5.11|) both the dependence on the 
Casimirs N and Cf as well as their proportionality tomp is in accordance to one's expecta- 
tions from the form of the corresponding HTL-resummed one-loop diagrams of Figs. [| and 
|6], respectively. 

However, in complete analogy to the scalar toy model of Sect. [II D|, we find that the 



magnitude of the corrections to the asymptotic masses are such that m 2 ^ + Sm 2 ^ drops 
to negative values for g > 1, which would give rise to tachyonic singularities in the semi- 



the naive strictly perturbative mass is 
For slightly higher values of rhn/T ~ g, 



perturbative entropy result (for N = 3 and \i 
again given by the shorter-dashed line in Fig. 
the same phenomenon occurs with + SM^. 

In the scalar model we have seen that including the NLO correction to the thermal 
mass in the approximately self-consistent form ( |2.61|) gives instead a monotonously growing 
function in g and very good agreement with the exact result in the N — > oo limit even for 
large g. For QCD, we define in analogy to Eq. ( |2.61| ) the NLA asymptotic mass through the 
quadratic equation 



g 2 (N + N f /2)T 2 g 2 NT 
6 v/2tt 



(5.12) 
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FIG. 11. The NLA entropy obtained by including Jmoo according to Eq. ( |5.12| ) for hard mo- 
menta k > A = \/2irTm£)CA in the pure-glue entropy. The central line in the shaded areas 
corresponds to c\ = 1, the two differently shaded areas to the bands ca = | . . . 2 and c\ = ^ . . . 4, 
respectively. The dotted line corresponds to a simple scalar model with constant mass m mod- 
ified according to Eq. ( 5.13|) such that it also contains the perturbative pure-glue result up to 
and including order g 3 ; the latter is displayed by the dashed line that leaves the plot already at 
m D /T ps 0.785. 



and similarly for M^. 

In contrast to the scalar case, however, where the thermal mass and its NLO correction 
was momentum-independent and therefore applicable for all momenta, the results (|5.11 ) 
apply only at hard momenta. Indeed, NLO corrections to thermal masses in QCD as far 
as they have been calculated turn out to be rather different at soft momenta: In Ref. |BU . 
the NLO correction for the plasma frequency of pure-glue QCD in the long-wavelength limit 
has been calculated with the result Sm^/m^ ~ —0.18\/~Ng, which is only about a third 
of the relative (averaged) correction of m 2 ^. Moreover, the Debye mass turns out to even 

x g logic /g), with recent lattice 



receive positive corrections |37] 5m D /rh 



simulations |38| yielding a rather large constant c. 



For this reason we choose to leave the HTL results for the soft gluonic propagators and 
self-energies completely untouched, and we implement the NLO correction to the asymptotic 
mass by introducing a cutoff scale that separates hard from soft momenta at a scale A = 
v27rTm#CA which is proportional to the geometric mean of the hard Matsubara scale 2nT 
and a soft scale ca^d- For momenta k < A we keep the HTL approximation and for k > A 
we take the thermal gluons to have the constant asymptotic mass fh 2 ^ of Eq. (|5.12| ). This 
completes the definition of our present next-to- leading approximation to the entropy Snla- 

In Fig. |ll|, the numerical result for pure-glue QCD with ca = 1 is given by the full line. 
The effect of varying ca in the range | . . . 2, which keeps A well in between the interval 
(m^, 2ttT) for all values of rhc/T, is displayed by the dark-grey band; the more extreme 
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variation ca = 7 • ■ ■ 4 extends the latter by the light-grey areas. We shall see, however, that 
in the eventual applications to QCD at temperatures a few times the transition temperature 
the resulting increase in our "theoretical error" will be still moderate when compared to the 
renormalization scheme dependences. 



For the sake of comparison with a simpler quasiparticle model, the dotted line in Fig. |TT 
shows the entropy of two interaction-free scalar bosons with constant mass 



*» = I 1 - *2^-L* I (5.13) 



7T 



which mimicks the NLA result ( |5.12 ), but adjusted such as to reproduce the perturbative 
QCD result up to and including order g 3 in this simpler model. 

In the fermionic quantities Sf and Af, which because of the absence of Bose enhancement 



are less sensitive to the soft scale, we implement the fermionic analog of Eq. (|5.12| ) by 
rescaling M 2 at all momenta for simplicity (as in the scalar case in Sect. II. D). 



C. Renormalization-group improvement 

In the HTL/HDL approximation, all the gluonic and fermionic contributions above de- 
pend on the numerical value of the HTL/HDL masses rh 2 D and M 2 , respectively, which are 
proportional to a s = g 2 /^. The latter is a renormalization scheme and renormalization 
scale (p) dependent quantity, and so are therefore our results for entropy and density. Fol- 
lowing Ref. Q we adopt modified minimal subtraction and assume that an optimal choice 
of the renormalization scale should be found around the scale of the Matsubara frequencies, 
27rT, or in the case of zero temperature and finite density around the scale of the diameter 
of the Fermi sphere, 2/i. After all, the hard thermal/dense loops are generated by hard ex- 
citations, as are in fact the NLO contributions to the asymptotic masses within HTL/HDL 
perturbation theory. 

Exactly as done in Ref. |TT| in a direct HTL resummation of the thermodynamic pressure, 
we put in the running of the coupling by hand and choose it to be determined by the 2-loop 
renormalization group equation according to 

^'m^V — imar) (5 ' 14) 

with L{jj) = log(/i 2 /A|jg) and 

A) = (lljV-2JV»/3, fa = (34N 2 - 13NN f + 3N f /N)/6. (5.15) 



1. Entropy 

At least at zero density, lattice results relate the QCD scale parameter Aj^g to the critical 
temperature T c , which in accordance with Ref. we choose as T c = 1.14Aj^g, both for 
pure-glue QCD and also for Nf 7^ 0, since lattice data indicate only a weak dependence of 
the ratio T c /Ay^ on the number of quark flavor. 
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FIG. 12. Comparison of the HTL entropy (full lines), the NLA results for ca = 
(dash-dotted lines) as well as ca = \ ...4 (gray dash-dotted lines), and the free entropy of bosons 
with mass (|5.13| ) such as to reproduce the correct perturbative plasmon effect (dotted lines), all 
with MS renormalization scale p, = ttT . . . 4irT, with the lattice result of Ref. [^0| for pure SU(3) 
gauge theory (dark-gray band). 



Putting ft = c^nT in Eq. ( |5.14| ) and assuming ~ 1 prescribes reasonably small values 
for a s and thus for rhu/(2iTT) and M/{tiT) for all T > T c so as to make it interesting 
to compare the above HTL and NLA expressions with nonperturbative results from lattice 
gauge theory. Indeed, we have found that, for <C 2nT and M <C ttT, the deviation 
from the free Stefan-Boltzmann result is small enough to make a semi-perturbative picture 
minimally tenable, although it is clear that the physics of the phase transition itself is 
completely beyond reach. On the other hand, the strictly perturbative results up to and 
including the order g 3 are such that entropy and pressure would be much higher than 
their Stefan-Boltzmann values, indicating a complete loss of convergence of strict thermal 
perturbation theory. 

In order to have some indication of the theoretical uncertainty involved, we consider, 

1 2 

2 . . . z. 



again as done in Ref. a variation of the renormalization scale by a factor of c- u 



For purely gluonic QCD, the lattice results involve the least uncertainties. In Ref. P0| , the 
thermodynamic potentials of pure SU(3) gauge theory have been calculated from plaquette 
action densities on lattices up to 8 x 32 3 for temperatures up to about 4.5T C and extrapolating 
to the continuum limit by comparing different lattice sizes. The lattice result for the entropy 
density is rendered in Fig. [12] by a grey band whose thickness is meant to give a rough idea 
of the errors reported in Ref. 
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Our result for the HTL entropy as displayed in Fig. [7| translates into a range of values 
bounded by the choices = irT (lower full line) and fl = 4nT (upper full line). This 
already gives a remarkably good approximation of the lattice result for T ^ 2T C , somewhat 
underestimating the values at higher temperatures. In all of this the parameter rho/T takes 
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on values in the range ~ 1 ... 2, where we have found the entropy of simple massive bosons 
to give only slightly larger results (cf. Fig. |?]). 

Now the HTL entropy contains only part (here 1/4) of the plasmon effect. The latter 
is completely included in the NLA entropy as defined after Eq. ( |5.12j ). In Fig. [12], Snla 
is represented by the area bounded by the black dash-dotted lines, where the lower one 
corresponds to the choice /i = 7rT and ca = 2, and the higher one to ft = 4nT and ca = 1/2. 
For this range of c\ the scale A remains well separated from both rhp and 2nT for all T > T c . 
Extending the range of Ca to Ca = 1/4 . . . 4 gives the area bounded by the gray dash-dotted 
lines. Although A now varies all the way from to 2nT, the error band is only moderately 
enlarged.^ 

Evidently, the NLA estimates based upon Eq. djj) do not move away too much from the 
HTL results, which first of all is what is required to make our semi-perturbative procedure 
tenable. What is more, the results show a surprisingly good agreement with the lattice 
results for temperatures greater than 2. . . 3 times the critical temperature. 

Recently, in Ref. ||1] the results of Ref. [40| have been reproduced within errors by using 
a renormalizat ion-group improved lattice action for temperature up to 3.5T C . The results 



of Ref. Ell] for the pressure are systematically higher by about 5. . . 2% for temperatures 



2. . . 3.5T C . For the entropy, which has not been extracted explicitly in Ref. [4~I], this would 
imply a result that is centered around the upper boundary of the grey band in Fig. 12 
for T « 3T C , and slightly flatter around T ~ 2T C , all with slightly reduced error bars. If 
anything, the agreement with our HTL and NLA results appears to be even a bit improved. 

Comparing finally with the entropy of free massive bosons with mass according to 
Eq. ( |5.13| ) such as to reproduce the correct perturbative plasmon effect, this is included 
in Fig. [12| as the band bounded by the dotted lines corresponding to Jjl — tcT . . . AttT. Since 
the renormalization scale dependence decreases with decreasing deviation from the Stefan- 
Boltzmann value, this band is rather narrow. It is also clearly in lesser agreement with the 
lattice data, which thus favor the momentum-dependent inclusion of NLO corrections to the 
thermal masses that follows from NLO perturbation theory and that we have modelled in 
our NLA estimates. 

In Fig. [13|, the central results for the HTL and NLA entropy (/2 = 2ttT and ca = 1) are 
displayed together with the results for Nf = 2 and 3. Only a rather weak dependence on Nf 
is found in this (greatly magnified) plot where the entropy is normalized to the free value, 
and T to the respective (^-dependent) T c oc Aj^g"- 

These results are in good agreement with recent lattice results [42] for Nf = 2 and their 
estimated extrapolation to the continuum limit and to the limit of massless quarks as we 



have already noted in Ref. [0 . In Fig. [13], a conversion of the lattice result to the entropy 
is included as a gray band, and, indeed, for T/T c ^2.5 our NLA estimate turns out to lie 
close to the center of the estimated error band of the lattice result. 



17 Although in Fig. 11 there was a noticeable increase in the error band for the NLA results when 
increasing the range of ca, this does not affect so much the total error because the lower bound, 
which corresponds to higher values of itid/T, is moved further down only by increasing ca, where 
the addition in the error band is small; the upper bound on the other hand corresponds to smaller 
values of rho/T, where the upward increase in the error is correspondingly smaller. 
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FIG. 13. Comparison of the HTL entropy (full lines) and the NLA result with c\ = 1 
(dash-dotted lines) for flavor numbers Nf = 0, 2, 3, all with the central choice of MS renormaliza- 
tion scale p, = 2ttT. The estimate of a continuum extrapolation of the lattice result for Nf = 2 as 



reported in Ref. 42.43] and its estimated error, converted to S/Ssb, is given by the gray band. 
Notice the blown-up scale of the ordinate compared to Fig. O . 



The further result that for Nf = 3 our NLA estimate for S/Ssb as a function of T/T c 
is approximately the same fits nicely to the recent lattice data for Nf = 3 |43|], which 
are consistent with a coincidence of the asymptotic values of P/Psb and also for S/Ssb- 
A more detailed comparison of our results with the lattice data, in particular at smaller 
temperatures, is hardly worthwhile in view of the large uncertainties associated with the 
extrapolation to the massless continuum limit .Q 



In our previous works |T5| , p!6| we have been considering a simple Pade-improved inclusion 
of the NLO asymptotic mass correction in place of the NLA form fl5.12| ). A comparison of the 
respective results shows that our estimate of the effects of an approximately self-consistent 
treatment of NLO corrections to the self-energies is fairly robust, with the main uncertainties 
coming from the choice of the renormalization scale. 



In a recent paper the authors of Ref. [44] have reported an extremely good fit of the entire 
lattice data using only the perturbative first-order correction to the pressure, a bag constant and 
a numerically integrated 2-loop /3-function. However, this agreement has been achieved with the 
lattice results which still contain finite-cut-off effects. In Ref. 42 ,43] , the size of the estimated 
correction for the continuum limit is given as +15 ± 5%. These corrections are essential for the 
good agreement with our results as shown in Fig. |l3| Conversely, the results of Ref. |44] remain 
even below the plot area of Fig. 13. 
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2. Density 



For nonvanishing chemical potential, where lattice data are missing to determine precisely 
the critical temperature or density in terms of Aj^g, we can nevertheless translate our results 
into functions of T/Aj^g and fi/Ay^ provided we choose the renormalization scale ft as 
a suitable combination of T and \x. If, as we have assumed, the spacing of the Matsubara 
frequencies, 2ttT, gives a good choice for the renormalization scale ft at zero density, it seems 
plausible to adopt the diameter of the Fermi sphere, 2/z, in the case of zero temperature. 
This choice of a relative factor of tt is particularly natural when considering the form of the 
leading-order result for the fermionic thermal masses, Eq. ( }4.21|) , where T and fi/n appear 
on equal footing. 

In Fig. n we give the numerical results for the quark density Mhdl at T = for N = 3 
and Nf = 3 as a function of fi/Aj^r for the range ft = /i . . . 4/z. In this case we do not attempt 
to include NLO corrections, for they do not contribute terms of order g 3 . NLO corrections to 
the hard fermion self-energy are in fact responsible for completing the plasmon effect at order 
g 4 log(g), but a complete calculation of the former would be needed to determine that part 
of the constant under the logarithm that comes from the spectral properties of quasiparticles 
rather than explicit order-g 4 interactions, which are dropped in the approximation J\f' = 0. 

The dashed line in Fig. [L4| gives the strictly perturbative result at order g 2 . The result 
corresponding to a simple quasiparticle model with mass M = is not included; from 
Fig. [10] it is clear that it is in between the HDL result and the order-g 2 one, and somewhat 
closer to the latter. 

The perturbative result up to and including order g A has been calculated by Freedman 
and McLerran []34| and by Baluni []45" |. However, it has been obtained in a particular 
momentum subtraction scheme. In order to convert this to the gauge-independent MS 
scheme, one should replace the scale parameter /z in Ref. |34| (M in Refs. |45|,[I|) according 
to 

// = ftexp {[(151 + 36a + 9a 2 )N - 40JV/]/[24(lliV - 2N f )]} , (5.16) 

where a is the gauge parameter used in the momentum subtraction scheme calculation. In 
particular for N = 3 and uniform chemical potentials one fmds[^] 



P/Po 



M/M = 1 



7i 



10.347- 0.536iV7 + N f log 



N f a s (fi) 
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+ (ll_|jV»lo g £ 
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7[ 



7.597- 0.369iV7 + N f log 

2 



N f a s (jj) 



TT 



+ 11 



TT 



(5.17) 



(5.18) 



19 The numerical coefficients have been assembled from Eq. (4.46) of Ref. |^5| using Eq. ( |5.16D , 
thus avoiding some unnecessary accumulated rounding errors that are present in the final results 
(5.14) and (5.15) of Ref. [45]. The actual error in the above numerical coefficients is probably 
about 1 in the next-to-last digit. 
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FIG. 14. The result for the quark density for Nf = 3 in the HDL approximation for p, = \i . . . 4/i 
(full lines) compared with the perturbative results at order g 2 (dashed lines) and order g 4 (dotted 
lines). 



With Nf = 3, this is included in Fig. 14" by the dotted lines^]. Although the perturba- 
tive order-g 4 result constitutes a substantial correction of the order-g 2 result, perturbation 
theory at zero temperature and high densities is clearly much better behaved than at high 
temperatures — the interaction terms are increased by less than 50% for /i ^ 2Aj^g (/x ^ 3A^g 
in the case of P/ Pq). 

On the other hand, the nonperturbative Mhdl result is rather close to the perturbative 
order-g 2 result, showing even a slight decrease of the interaction contribution compared to 
the latter. The Mhdl result contains already a fraction of the coefficient of the a 2 log(a s ) 
term, together with a subset of the true higher-order contributions. It would be interesting 
to see how an NLA calculation, which would complete the g 4 log(g) coefficient, compares 
with the perturbative order-g 4 result. We intend to investigate that in future work. 



VI. CONCLUSIONS AND OUTLOOK 

We have shown that it is possible to perform a resummation of HTL's which is free 
of overcounting and UV problems through approximately self-consistent calculations of the 
thermodynamical functions of QCD, without the need for thermal counterterms. The two- 
loop skeleton approximation for the free energy reduces to effectively one-loop expressions 
for the entropy and the density but with dressed propagators. With the latter approximated 



20 In Fig. 4 of our previous publication Ref. [16|, the perturbative order-g 4 result was not correctly 
included because of an incomplete scheme conversion. 
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by the HTL /HDL propagators we reproduce correctly the leading-order interaction termsQ. 
In fact, the latter can be expressed entirely in terms of the asymptotic thermal mass. The 
so-called plasmon-effect contributions of order g 3 on the other hand are only partly ac- 
counted for by HTL self-energies and propagators; the remaining contributions arise, rather 
unconventionally, from NLO corrections to the self-energy of hard particles at the light-cone 
as given by standard HTL perturbation theory. 

This is to be contrasted with a direct HTL resummation of the one-loop pressure |TT 



There the plasmon effect is contained completely in the soft contributions, whereas the 
leading-order interaction terms are over-included and corrected only in a (very complicated) 
two-loop calculation. 

We would like to recall that while the HTL and NLA approximations that we have con- 
sidered are manifestly gauge independent, this is not the case for our starting point, the 
self-consistent ^-derivable two-loop order approximation itself. The corresponding gap 
equation would involve unphysical gauge dependent features as well as an incomplete low- 
est order /3-function both of which enter at order g 4 , i.e. beyond the perturbative accuracy 
of a two-loop approximation. These are automatically dropped in our current approxima- 
tions. Further improvements, beyond our HTL and NLA approximations would require to 
also improve upon the self-consistent two-loop approximation. In order to achieve gauge 
independence in (approximately) self-consistent resummations one should obviously turn to 
approximations which include dressed vertices, using for instance the formalisms which have 
been developed long ago by de Dominicis and Martin [17] and also Freedman and McLerran 
[ 34fl . The strategy, in principle, would be to include vertex corrections, together with in- 
creasingly better approximations for the quasiparticle propagators. That is, with increasing 
number of loops in $, the building blocks in this scheme — the self-consistent propagators and 
vertices — should be also improved. However, a practical implementation of such a scheme 
in the case of nonabelian gauge theories seems to be hopelessly complicated. It is therefore 
gratifying that the approximate propagator renormalization that we have presented turns 
out to be already a good approximation. 

In the expressions that we use for the entropy and density the main contribution comes 
from the vicinity of the light-cone where hard thermal loops remain accurate also at hard 
momenta and provide the asymptotic masses. We have proposed a procedure of including 
NLO corrections through approximately self-consistent corrections to the thermal masses of 
the hard excitations only. The NLO corrections to the asymptotic masses can be calculated 
more accurately by means of standard HTL perturbation theory, the details of which are 
postponed to a forthcoming publication. 

The numerical evaluation of our results combined with a two-loop renormalization group 
improvement turn out to compare remarkably well with available lattice data at zero quark 
chemical potential, which supports the picture according to which much of the effects of 
the interactions in the quark-gluon plasma can be adequately described by means of weakly 



an attempt has been made to resum the HTL self-energies directly on the level of the 
skeleton representation of the free energy. However this relies on an arbitrary modification of the 
functional <1> which, although it yields the correct g 2 effects (by construction), does not respect the 
correct combinatorial factors and thus violates the proper counting of the higher-order diagrams. 



21 InRef. [46 
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interacting gluonic and fermionic (HTL) quasiparticles. 

Extensions of the present work which are in progress concern the evaluation for general 
fi > and T > 0, and the integration of entropy and density to the thermodynamic pressure 
P(/i, T), similarly to what has been done in simple quasiparticle models in Ref. j|7|] (see also 
fl6|j). Maxwell's relations, which constitute the corresponding integrability conditions, are 
satisfied up to and including order g 3 upon inclusion of the NLO contributions; beyond that 
order they give constraints on a possible renormalization-group improvement and it seems 
interesting to further pursue the present approach of combining the physical content of the 
perturbatively derivable hard thermal/dense self-energies with nonperturbative expressions 
for entropy and density, which in self-consistent two-loop order approximations only depend 
on the spectral properties of quasiparticle excitations. 
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APPENDIX A: THE PLASMON EFFECT IN THE QCD ENTROPY 

In this Appendix, we shall explicitly verify that our approximation for the entropy (cf. 
Eqs. (|3.10|) and (|4.10| )) contains indeed the right perturbative correction of order g 3 . Recall 
that £3 involves two types of contributions: the LO entropy of the soft gluons (longitudinal 
and transverse; cf. Eqs. ( |3.31|) and (|3.32| )), and the NLO entropy of the hard particles 



(transverse gluons and fermions), as determined by the corresponding NLO self-energies on 
the light-cone (cf. Eqs. (|3~30D and flL2p). 

Our strategy will be as follows: In Sect. A.l, we shall rewrite the soft gluon entropy in a 
way which will be convenient later. Then, in Sects. A. 2 and A. 3 we shall compute the NLO 
self-energy SUt of a hard transverse gluon, and the corresponding contribution 5St to the 
entropy. This will complete the derivation of the plasmon effect for a purely gluonic plasma. 
The extension to a plasma with fermions will be finally considered, in Sect. A. 4. 

1. The entropy of soft gluons 

From Eqs. ( |3.31 ) and ( 3.32 ), the order-g 3 contribution of the soft gluons reads 



S s 3 o(t = - \ { Im [log(l + A>n) - UDo] - ImflRe (D - D )} 



S ( 3 a) +AS 3 , (Al) 



where = {dP z /dT)\^ D (cf. Eq. flO§)), and AS 3 is defined in Eq. ([T43]). In Sect. 



Il l D3| , we have mentioned that AS3 has been numerically found to vanish, because of a 
compensation between the electric and the magnetic contributions to Eq. ( |3.43| ). In what 
follows, however, we shall not use this information, but rather consider separately these 
electric and magnetic contributions, and show how they combine with the corresponding 
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contributions to the NLO entropy of the hard particles 5S. Specifically, we shall verify that 
the identity in Eq. ( 3.36 ) holds separately in the electric and the magnetic sector. 

To this aim, it is convenient to rewrite Eq. ( |3.43| ) in a slightly different form by using 
ImllReZ) = Im (III)) — RelllmZ), and then integrating the first term: 

d 3 k r duo T Ai, ,a -|(o)\\ _ ™b f d 3 k (\ 1 \ rh% 



(2,)3 J 2™ Im ^ ~ B ™» = 2 /(Sr)» (*3~ j = 8? ' (A2) 

This yields (with p L>T = 2Im J D LiT , cf. Eq. Q): 

A5 3 = i {/3L(ReriL - m^) - 2(p T - Refl r } = AS< 3 > + a4 3) , (A3) 

where we have also used the following "sum- rule" (cf. Eq. ( |3.6| )): 

du pl{oj, k) 1 1 



2n uj k 2 k 2 + m 2 T 



D 



(A4) 



Given the complicated structure of the HTL self-energies and spectral functions, the integrals 
in Eq. ( |A3| ) cannot be further evaluated in closed form. But this is actually not needed: 
indeed, the cumbersome terms in these expressions will be shortly shown to cancel against 
similar terms in £3 rd , the order-g 3 contribution of the hard particles, to be computed below. 



2. The NLO gluon self-energy 

We shall now compute the NLO self-energy contribution SUt of a hard transverse gluon. 
This is determined by the effective one-loop diagrams in Fig. ^ where one of the internal 
lines is a soft gluon (L or T) with the HTL-dressed propagator D — D (the subtraction of 
the free propagator D is ensures that the loop integral is saturated by soft momenta). The 
other line in each of these diagrams is hard and transverse, and therefore undressed. 

We are interested only in the transverse projection of SU^: 

5U(p) = 5U T (p) = ~ - m flly(p). (A5) 

We write 5Ut = 5U l + 511', where the upper indices refer to the soft internal lines in 
these diagrams, and compute only the longitudinal contribution 511' in more detail. (The 
calculation of the transverse contribution is completely analogous.) This involves two of the 
diagrams in Figs. |j: the tadpole 5U l a and the non-local diagram 5111, . The tadpole gives 

Sn' a = -g 2 Nj[dk](D L (k)-D^(k)) 

where the Matsubara sum in the first line has been performed by using the spectral repre- 
sentation ( |3.6|) , and in the second line we have replaced n{k Q ) ~ T/k (as appropriate at 



d 4 k 

W) 
g 2 NTm D 

Att 



r d k 

~9 2 N Jj^^p L {ko,k)n(k ) 

(A6) 



56 



soft energies), and then performed the energy integral with the help of the sum-rule ( |A4| ). 
The final result in Eq. (A6) is indeed of order g 2 Trfi£, ~ g 3 T 2 , as expected. 
The non-local diagram in Fig. ^.b yields: 

(6n l b )..(p) = -2^f J[dk](2 Po + k f Dlf(p + k)(D L {k) - 4 0) (A;)), (A7) 

where D^(q) is the free magnetic propagator in the Coulomb gauge, 

d$\q) = (A8) 

yo — h 

and the factor 2 in front of the integral reflects the two possible ways to choose the soft 
longitudinal line among the two internal lines in the original one-loop diagram. 

The transverse projection of Eq. ( |AT| ) involves \ —p l p')(Sij — qify), where q = k + p. 
Since p ~ T, we have qi = j^^i — V%i while the integral in Eq. ( |A7j ) will be eventually 
dominated by soft k momenta. In what follows, we shall often perform such kinematical 
simplifications relying on the fact that k <C p. With this simplification, the product of the 
transverse projectors above reduces to the identity, so that 

5U l b {p) = ~g 2 N J[dk}(2p + k ) 2 D (p + k)(b L {k) - D^\k)). (A9) 

To perform the Matsubara sum over ko we need the following sums (with q = k + p): 

r£B„( P + k)(D L{k ) - Df\k)) = fp M k) pM ", fa) "f o) , (Aio) 

V J 2n J 2n kn — On + Vn 



k -q + p 
fdk [dq n(q ) - n(k ) 



Tj2k D {p + k)(b L (k)-DP(k)) = [p- f^k p L (k) Po (q) 

ko v J 2tx J Z7r k - g + Po 

T^k (k +p )D {p+k)(b L (k)-D^(k)) = f d ^k q p L (k) P0 (q) f 0) - n } k0) . 
This finally yields, for the retarded self-energy, 

H(p) = -g*N f^p L (k) P o(q) M + Sp k + k ogo ] , n(9o) 7 n( y. . (All) 

J (Zn J 4 J In k — q + p + ie 

To compute the entropy ( |3.30| ) we need the light-cone projection of the real part of this 
self-energy, 

Re5l4(p = p). Note that in previous calculations of the damping rate, it was 
rather the imaginary part of this same self-energy which was required The calcula- 

tion of the imaginary part is easier since the LO contribution ~ g 2 T 2 can be immediately 
extracted from Eq. ([A lip by neglecting n(q ) ~ 1 against n(k ) ~ T/k Q ^> 1, and keeping 
only the large external momentum 4p^ in the numerator. This, together with 

Im — = -ir5(k - q +po) ~ -n5(k - k cos 9 + po - p), (A12) 

h-qo + Po + ie 



leads to the following, standard, result for the longitudinal part of the damping rate [48 
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_ lm5U[(p =p) g 2 NT r d 3 k rduj 
li = o = — o — TTT^ / — Pl\u,k)o{w- kcosV). (A13) 



2p 2 J (2vr) 3 J 00 

If we perform, however, the same simplifications on the real part, then the would-be LO 
result turns out to vanish, by parity (with P denoting the principal value): 

(Re5Ui(p =p)) . = g 2 NT f^L f—p L ^k)P ]— = . (A14) 

v 'naive J {ZlTj 6 J TTUJ UJ — k COS V 

In fact, this is only to be expected: the terms in Eq. (|A14 ) are formally of order g 2 T 2 , while 
we know that Re 511 should be rather of order g 3 T 2 . Thus, in order to extract the leading 
contribution to Re 511 from Eq. ( All ) one has to push the kinematical approximations one 
step further as compared to the damping rate. In particular, we need the expansion of the 
statistical factors in Eq. (|A11| ) to LO and NLO order: 

n ^ r \ T 2n (<?o) + 1 / Aik n 

n{k ) - n(q ) ~ . A15 

k 2 

We shall denote by Re 511^ the contribution coming from T/ko, and by Re<5Il| 2 the re- 
maining one due to (2n(q ) + l)/2. These quantities will be evaluated at po = p, so they are 
functions of the three-momentum p alone. We have: 

Re SU iM = fNT /£. Mk) f% [V + + M„] P ^ 

d 3 k r dk, 



Qo + P 



= g2NT Ij^j2^ h{k){2p + ko)2ReD ° {P + kh (A16) 

where in the second line we have identified the (retarded) free propagator via its spectral 
representation. In Re^II^, we can restrict ourselves to the LO term 4p 2 in the denomi- 
nator, and to the positive-energy pole go = |p + k| ~ p + kcos9 in the spectral function 
Po(Qo, |p + k|). This yields: 

Re H M , V WO + 1) • (AIT) 



3. The NLO entropy of hard gluons 



Let us consider first a purely gluonic plasma, in which case the hard gluon self-energy 
SUt = 5II Z + 5IP, is all what we need to compute the NLO entropy <5>3 ard = SS l + 5S l . As 
before, we focus on the longitudinal contribution SS l ; by inserting Eqs. (|A6|) , (|A16|) and 
( |ST7D in Eq. flOOP , we obtain 5S l = 5S[ + 5S l 2 , where: 



SS[ = —N, 



d A p 
—g 2 NN g T 



Po(p) 
d A k 

p L (k) 



dn(p ) 



dT 
PL{k) 



su l n + Re mi 
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(2tt) 4 2k 



k 
d_ 
dT 



d A p dn(p 



(2tt) 4 dT 
Refl L {k ,k). 



p (p) {2p + k ) 2 ReD (p + k)-l 



(A18) 
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In writing the last line above, we have identified the one-loop contribution to the self-energy 
of the soft longitudinal gluon due to hard transverse gluons. To the order of interest, this is 
precisely the HTL 11^. The second piece 5S l 2 of the entropy reads: 



SSl 



~N 9 
9 2 NN g 

2 dT 



d 3 p 1 dn(p) 



(2tt) 3 p dT 
d 3 p 



Re5U l b2 (p) 



(2tt] 



3 (2n(p) + 1 



dn(p) f d 4 k PL(k Q ,k) 



d A k 



dT 

PUh, k) 



(2 



7T 



k cos 9 



(2vr) 4 k -kcosO 



-N, 



(Tml) 



J 2vr 4 2k n art V D Ji 



(2tt) 4 2k ofvy-^-'^jr ( A19 ) 

In going from the second to the third line above, the following chain of identities has been 
used (see also Eq. ( |3.16| )): 

d 3 p 



2g z N 



2g z N 







dT J (2tt 



y dT J (2tt) 3 dp dT^ D > 



3 n(p)[l + n(p)\ 
dn 



d 



dp 



(A20) 



Furthermore, in writing the last line in Eq. (|A19 ), we have identified the longitudinal HTL 
RelL; as follows (compare to Eq. (|3.14|) ): 



Refl L (uj,k) 



-m 



D 



dfl k cos 9 
An uj — k cos 9 



(A21) 



By adding Eqs. ( |A18|) and (|A19 ), we finally deduce the following expression for the 
longitudinal piece of the NLO entropy: 



5S l 



N 9 T^ 

9 g T 



d^k Phik) 



(2vr) 4 2k Q 



N n 



d^k pi(k) 



(2vr) 4 2k 



(Re IT, 



m 



D 



dm 2 D NgTriD 
dT 8tt 



- AS?, 



(A22) 



with AS^ as defined in Eq. ( |A3|) . An entirely similar calculation shows that the remaining, 
transverse, piece <5iS* cancels against the transverse contribution AS^ to A1S3, Eq. flA3|) : 



<5<S* + AS 



(3) 



0. 



(A23) 



That is, the total contribution of the soft transverse gluons to the plasmon effect cancels 
away, as it should. 

All together, Eqs. ([All), (A3), ( A22 ) and ( A23 ) provide the expected result for the order- 
g 3 effect in the entropy of the purely gluonic plasma where T(dTrh 2 D ) 

-1 3 ^-iivC-i 2 



2m 2 D : 



csoft 
^3 



+ s 



hard 



N g m 6 D T dm 2 D N g rh D 



12tt 



dT 8 



7T 



N g rh 3 D 
3tt 



(A24) 



Moreover, it can be easily recognized that Eqs. ( A22|) and (|A23|) are equivalent to the 
longitudinal, and, respectively, transverse components of Eq. (|3.36| ), as they should. 
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4. Adding the fermions 



The previous results are easily extended to a QCD plasma with fermions. The entropy 
^hard j n ^ g cage i nvo i ves a i so fa e NLO self-energies of the hard fermions, <5E± : 

d 4 P ( AT dn( Po ) , _ NTJ _ , nATAT df 



5 3 hard = -/^|iV 9 ^^po(p)Rem T + 2NN f ^J2 Ps (p)Re5Z s y (A25) 

Once again, we focus on the contribution 5S l of the soft longitudinal gluons, and use the 
integral over the hard momentum p in Eq. QA25 ) to reconstruct the HTL 11^. Here, this 
involves both a hard gluon loop and a hard fermion loop, which enters via the self-energies 

For instance, the fermionic analog of R,e5U l b2 , Eq. ( |A17|) , reads: 

Re5E±2(p) = -^ [1 - 2/±(p)] i (2^ fc -A:cosr (A26) 

which, when inserted into Eq. ( |A25| ), determines the following contribution to the NLO 
entropy (compare to Eq. (|A19|) ): 

5S > = 9 N ° Jj2^i N{1 + 2n) df + T S (1 ~ 2fs) df\J k -kcose 

~ 2 9TV im Wi(2vr)^o-A;cosr (A27) 

(We have used here (1 - 2f)d T f = d T [f(l - /)] = -d T (Td k f), together with Eq. Q for 
the Debye mass.) This is formally the same result as for pure glue, Eq. ( |A19|) , except that, 
here, rho is the full HTL Debye mass, which includes contributions from fermions. 

Similarly, the other contribution 5S[ preserves the form in Eq. (|A18|) where, however, it^ 
is now the full HTL in a theory with fermions. Thus the final result in Eq. (|A22|) is formally 
unchanged, but it now applies to a QCD plasma with fermions, for which T(dTTh 2 D ) = 2m\ 
(cf. Eq. (HQ)). 

Consider finally the order g 3 effect in the quark density: as explained in the main text, 
this comes entirely from the NLO corrections 5H± to the hard fermion self-energies, and, 
more precisely, from the longitudinal sector alone (the soft transverse effects eventually 
cancel, as in the case of the entropy). Thus, A/3 = 5Af l , with 8M l given by the same 
equations as above, except for the replacement of the temperature derivatives by derivatives 
with respect to p. Thus 5N l = 8N[ + 6N%, where (cf. Eqs. ( \KWj ) and (|Mg)): 

d 4 k f>L d 



W '- N ' T 1(2^^8-^^' < A28 » 



and 



x Kf i_ N 9 Td ™ 2 D fd 4 k p L (k ,k) _ NT [d 4 k p L d , . 2 x 
5N * ~ -2~diT J (2k)* ko-k cosO ~ 9 7 (2^)4 2k~ ^i ReIlL " mD ) ■ (A29) 

As in the entropy, the non-local terms involving Ren^(fc) cancel in the sum of the two 
contributions above, and we are left with the following simple expression: 

dml r d A k p L {k) N g Tmlm D 

Ms ~ N9T ^p-Jj2W^kV ~ Airp ■ (A30) 
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